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We describe the gas of ultracold fermionic atoms by a functional integral for atom and molecule 
fields. The crossover from Bose-Einstein condensation (BEC) to BCS-type superfluidity shows uni- 
versal features in terms of a concentration parameter for the ratio between scattering length and 
average interatomic distance. We discuss the relevance of the Yukawa coupling between atoms 
and molecules, establish an exact narrow resonance limit and show that renormalized quantities 
are independent of the Yukawa coupling for the broad resonance, BCS and BEC limits. Within 
our functional integral formalism we compute the atom scattering in vacuum and the molecular 
binding energy. This connects the universal concentration parameter to the magnetic field of a 
given experiment. Beyond mean field theory we include the fluctuations of the molecule field and 
the renormalization effects for the atom-molecule coupling. We find excellent agreement with the 
observed fraction of bare molecules in 6 Li and qualitative agreement with the condensate fraction 
in 6 Li and 40 K. In addition to the phase diagram and condensate fraction we compute the correla- 
tion length for molecules, the in-medium scattering length for molecules and atoms and the sound 
velocity. 
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I. INTRODUCTION 

Ultracold fermionic atoms are developing into exciting 
experimental laboratories for the understanding of com- 
plex many body quantum physics. Collective phenom- 
ena as the phase transition to low temperature super- 
fluidity can be investigated with controlled microscopic 
physics. Recent experimental progress Q, 0, 0, 0, HJ in 
the crossover region between a Bose-Einstein condensate 
(BEC) 0,0 of molecules and the condensation of corre- 
lated atom pairs similar to BCS-superconductivity 0, EJ 
reveals the universality of the condensation phenomenon, 
as anticipated theoretically [Tol HI 111 111 111 ITU . 

At first sight the understanding of the crossover region 
in presence of a Feshbach resonance seems to depend on 
complex atomic and molecular physics and therefore on 
many parameters. The crucial role of molecular bound 
states involves the detailed knowledge of binding or reso- 
nance energies as well as atomic scattering amplitudes in 
presence of a homogeneous magnetic field B. In particu- 
lar, one needs the effective coupling between the bound 
state and the atoms, usually encoded in the matrix ele- 
ment of the Hamiltonian between "bound and open chan- 
nels" in the space of two-atom states. The many body 
physics depends in addition on the temperature and den- 
sity of the thermodynamic equilibrium state. Further- 
more, the experimental observations are influenced by 
the geometry of the trap which is used for the confine- 
ment of the atom gas. 

In this paper we underline the universal aspects of the 
equilibrium state of ultracold fermionic atoms. In the 
limit of a homogeneous situation we argue that "macro- 
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scopic" observables can be expressed in terms of only 
three dimensionless parameters: 

(1) The "concentration" c describes the ratio between 
the in-medium scattering length and the average distance 
between two unbound atoms or molecules. For small neg- 
ative c the molecular resonance becomes unimportant - 
the gas is approximated by fermionic atoms with a small 
pointlike interaction. The phase transition to superfluid- 
ity is of the BCS type and we will call this region in pa- 
rameter space (|c| « 1, c < 0) the "BCS regime". In con- 
trast, for small positive c the phase transition and the low 
temperature physics is well described by a gas of bosonic 
molecules. This limit will be denoted by "BEC regime" - 
it shows the phenomenon of Bose-Einstein condensation 
very similar to a gas of bosonic atoms. The perhaps most 
interesting region is the "crossover regime" between the 
BCS and BEC regimes for large c or | c 1 1 < 1. All three 
regimes are continuously connected as c _1 is varied from 
large negative to large positive values. For low enough 
temperature T a transition to a superfluid state with con- 
densation of atom pairs or molecules occurs for all c. The 
critical temperature T c for the phase transition to super- 
fluidity clearly shows (cf. fig. the universal character 
of the crossover. 

(2) The details of the crossover physics depends on 
a dimensionless Yukawa or Feshbach coupling which 
describes the interaction between atoms and molecules. 
This coupling depends on the density hj, cx nT 1 / 6 . In the 
"narrow resonance regime" of small the precise value 
of the coupling is not important for the location of the 
phase transition. A universal "narrow resonance limit" 

— ► exists for the "symmetric phase" for T > T c . For 
this limit we find a nontrivial exact solution of the many 
body problem which is described by a suitably adapted 
mean field theory. Small deviations from this universal 
limit lead for many (but not all) quantities to a con- 
trolled expansion for small h$. For the "broad resonance 
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regime" with large the precise value of does not 
affect the observables that are connected to renormal- 
ized fields and parameters. In fact, for large strong 
renormalization effects lead to a new form of universality 
for large couplings as observed previously in relativistic 
systems with a strong boson-fermion "Yukawa coupling" 
[lfij . Our general universality hypothesis 1 states states 
that the "macroscopic quantities" can be expressed in 
terms of only two couplings c and h$ . 

(3) The third parameter is finally the temperature in 
units of the Fermi energy, T = T/cf, £f = k F /2M, 
kp = (Z^n) 1 ^ (with M the mass of the atoms). Us- 
ing the parameters c, and T we obtain results for 
dimensionless ratios and couplings for arbitrary densities 
n. For a given c and T the only implicit dependence on 
n arises via the weak dependence oc n -1 / 6 and drops 
out in the limits of "enhanced universality" for broad and 
narrow resonances and the BEC or BEC regimes. 

No further details of the "microphysics" are needed for 
the macroscopic quantities. In this sense the description 
becomes universal. In the language of quantum field the- 
ory the parameters c and h$ describe relevant couplings 
for the long distance physics. 

At this point we may emphasize that the universal 
aspects discussed here are not related to the "universal 
critical behavior" near a second order phase transition. 
For example, the divergence of the correlation length for 
T — > T c is governed by a universal critical exponent - this 
has not yet been computed very accurately here. The dis- 
cussion of universality in the present paper rather con- 
centrates on the parameters which are of relevance and 
concerns "non-universal quantities" in the sense of criti- 
cal phenomena. 

Our results are based on a systematic formulation of 
the interacting many body system in terms of a func- 
tional integral. A detailed derivation and discussion of 
this formalism, together with an extension to the inho- 
mogeneous situation of a trap, can be found in ^^|. We 
carefully address the issue of renormalization or, more 
precisely, the relation between "microscopic parameters" 
and "macroscopic observables". Expressed in terms of 
our renormalized parameters c, and T all loop inte- 
grals are ultraviolet convergent such that the actual value 
of the physical microscopic cutoff scale plays no role. Also 
the mass of the atoms drops out. 

As one of the important advantages of our functional 
integral formulation we can compute both the "physics in 
the vacuum" and the atom gas at nonvanishing density 
and temperature in the same framework. The vacuum 



The universality here differs from the universality argument given 
by Ho Il7l . In particular, in the crossover region the interaction 
is not pointlike and does not become scale-free for the limit of 
diverging scattering length oc |c|. The Yukawa coupling intro- 
duces a further mass scale which manifests itself in the dimen- 
sionless parameter h^. Ho's argument can be applied, however, 
in the limit c~ 1 — > 0, h& — > oo. 



physics corresponds to appropriate low density limits and 
concerns properties of two-atom systems like the cross 
section and molecular binding energy. These quantities 
are directly accessible to experiment. Combining the ex- 
perimental information about the two-atom system with 
knowledge of the density (fcp) we can relate the param- 
eters c and htj, to measurable quantities. In particular, 
the concentration c is related to the magnetic field, or, 
more precisely, to its "detuning" AB = B — B from the 
vacuum Feshbach resonance at Bq. 

The functional integral treats the fermionic fluctua- 
tions of unbound atoms and the bosonic fluctuations of 
the molecule or di-atom field on equal footing. It is 
therefore an ideal starting point for the inclusion of the 
molecule fluctuations. These fluctuations are important 
for the quantitative understanding of the phase transition 
for a broad Feshbach resonance. For a large dimension- 
less Yukawa coupling h$ the renormalization effects are 
crucial. In fact, the importance of the molecule fluctua- 
tions is directly related to the size of - for h$ — > a 
nontrivial exact "narrow resonance limit" can be estab- 
lished where the molecule fluctuations can be neglected 
and mean field theory becomes valid. The present exper- 
iments in 6 Li and 40 K probe, however, broad resonances 
with large h<j, for which the inclusion of the bosonic fluc- 
tuations cannot be neglected. 

We present a first computation of the correlation 
length for the molecule field. It turns out to be larger 
than the average distance between two atoms for a rather 
broad range of temperature above T c , indicating the im- 
portance of collective effects. We also compute the effec- 
tive scattering length for molecule-molecule scattering. 
Again, it turns out to be large in the crossover region. 
Furthermore, we compute the sound velocity in the su- 
pcrfluid phase. 

Our functional integral approach permits a unified 
view of many previously obtained results. In fact, vari- 
ous other methods have already revealed many aspects of 
the results of our investigation. Strinati et. al. propose 
a purely fermionic description of the crossover problem 
based on a physically motivated choice of diagrams in the 
frame of the operator formalism 0, [23, 0, l22l I23I ] . As- 
suming that the chemical potential satisfies fi/T — > — oo 
in the BEC limit, they can show by explicit calculation 
[22l| the appearance of a "molecular" and a condensate 
particle density in the BEC limit. A similar approxima- 
tion scheme is advocated in [24| . 

Other approaches formulate the problem microscop- 
ically in terms of a Yukawa theory for a coupled sys- 
tem of fermionic unbound atoms and bosonic molecules 
IllplllllMllSllilsllllIp,!!! Illjp. In part, 
this is motivated by a functional integral [26l . l33l . l34j . Ap- 
proaches of this type are well suited for dealing with non- 
local interactions. Another route of achieving this task 
is followed in |37l l38| . A comparison of purely fermionic 
vs. bosonized approach has recently been performed |39| . 
Working with the microscopic parameters appropriate for 
40 K, almost no difference is found between the two ap- 
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proaches. This is consistent with our investigation, clas- 
sifying 40 K in the broad resonance regime where the two 
approaches become indeed equivalent. 

A systematic study relating on the two body level the 
parameters of the microscopic Yukawa theory to binding 
energies and two-atom cross sections in vacuum includes 
the UV-renormalization [33L 133. |3M H(j . The investiga- 
tions of Feshbach resonances (see |2a| for K) reveal a clear 
tendency that the Yukawa coupling extracted "naively" 
from scattering properties is much too large in order to 
account for the large condensate fraction ttc advocated 
by observation, in particular for 6 Li used by many experi- 
mental groups 0,0,E3)0)E3- This has triggered recent 
work on an effective Yukawa coupling and its momentum 
and energy dependence [3(| |43|, |44( • 

The plan of the paper is as follows. In sect. [H] we 
briefly discuss the formal ingredients of our functional 
integral formalism. It is based on the effective action to 
be computed for our two-channel model of the crossover 
system. A more extensive description of these issues can 
be found in ^^|. In sect. IIIII we state our notion of 
universality more precisely. The aspects we consider dif- 
fer from those analyzed by other authors, e.g. 0. We 
state universality here in a general field theoretical con- 
text, where for given external parameters n and T only 
two relevant renormalized couplings describe the system, 
namely c and h^. In sect. II VI we briefly address enhanced 
universality in various limits. These limits are discussed 
in more detail in sect. IVIIII when the necessary formal 
framework has been provided. 

In sect. 0we turn to the important difference between 
bare (microscopic) and dressed molecules (or "macro- 
scopic" di-atom bound states). We compute the renor- 
malization consant that relates the bare and dressed 
molecule fields. Sect. IVII turns to the effective potential 
for the molecule field. Here we state our approximation 
scheme beyond mean field theory. 

In sect. I VIII we compute the correlation length for the 
molecule fluctuations. It diverges at the critical temper- 
ature and remains large for a rather wide range, typically 
as long as T is smaller than 1.5T C . A large correlation 
length is a clear sign for collective effects. We also com- 
pute the in-medium scattering length for molecules and 
the sound velocity in the superfuid phase. 

In sects. HXIXI we relate the universal couplings c, h<p to 
the parameters of Feshbach resonances, in particular the 
detuning of the magnetic field B — Bq and the response 
of the level splitting to variations of the magnetic field, 
p,. Our functional integral formalism enables us to com- 
pute directly the observable quantities of the two-atom 
system, like the molecular binding energy cm or the scat- 
tering cross section between the unbound atoms. This is 
achieved in an appropriate limit where both the temper- 
ature T and the density n vanish as shown in sect. IIXI 
Sect. E] establishes a direct link between the universal 
concentration parameter c and the magnetic field B and 
determines the Yukawa coupling h<f,. 

In sect. IXII we turn to our results for the crossover 



phase diagram for particular Feshbach resonances ex- 
plored by recent experiments. We first address the frac- 
tion of bare molecules &b that has been measured by 
laser excitations to a higher level p|. We find excellent 
agreement with the observations in a range covering five 
orders of magnitude for fi^ on both sides of the Feshbach 
resonance (cf. fig. OJ- The condensate fraction fie relates 
to dressed molecules. Here we find qualitative agreement 
with the experiments in 40 K Q and 6 Li 0, whereas a 
quantitative comparison would need a more detailed re- 
lation between the experimental observable and ilp. In 
sect. IXIII we present our conclusions. Four appendices 
deal with the renormalization effects for the Yukawa cou- 
pling |Q[J, the detailed form of the gap equation J5J), the 
dispersion relation for the molecules in vacuum and 
the atom scattering in vacuum (|P|) . 



II. EFFECTIVE ACTION 

The ultracold gas of fermionic atoms in the vicinity of a 
Feshbach resonance can be idealized by two stable atomic 
states denoted by a two component spinor %p. (For the 
example of 6 Li these states may be associated with the 
two lowest hyperfine states |1), |2).) The molecular state 
responsible for the Feshbach resonance can be treated as 
a complex bosonic field </>. In our approximation this bo- 
son is stable for negative binding energy and can decay 
into a pair of fermionic atoms for positive binding en- 
ergy. Our functional integral formalism is based on the 
euclidean microscopic action 



f r A 
S R = Idx V f (9 T - — - v)ip 

^-(j>*i() T ei{}- ^etp" 



(1) 



This is a simple model for fcrmions with Yukawa cou- 
pling to a scalar field 4> and a pointlike "background " 
interaction (or "four-fermion vertex" ) A,/, . The parame- 
ter a plays the role of an effective chemical potential 2 . 
The different prefactor for the coupling between a and 
the atoms resp. molecules accounts for the fact that the 
molecules have atom number two. Similarly, their mass 
is twice the atom mass such that they have a nonrela- 
tivistic kinetic energy p 2 /AM. Further motivation and 
details of our functional integral approach can be found 
in |lSj . Our units are h = c = ks = 1. 

The bare parameters &a , and A^, have to be fixed by 
appropriate observable "renormalized" parameters. The 



2 We treat here only homogeneous situations such that the more 
complete formalism of lid is not needed. We also have not 
bosonized the four-fermion vertex ~ . In the language of |1S|I 



this corresponds to the formal limit m 
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quadratic term ~ <j)*4> involves the "bare" binding energy 
(Da) which typically depends on the magnetic field. The 
Yukawa coupling accounts for the coupling between 
the single atoms and molecules, e a p — — ep a , e\2 = 1. 
For h,f, — > the molecular states decouple, giving rise 
to a nontrivial exact "narrow resonance limit", cf. sect. 
IV111AI 

For a thermodynamic equilibrium situation the parti- 
tion function is represented as a functional integral with 
weight factor e" s , with S the Euclidean action 



, [U] = I XtyZty exp { - S B [V, a, $\ (2) 



dx\ 



Here Zb is defined as a functional of the source term 
j(f>(x) and we will set jj, = at the end of the computa- 
tions. 

Our Yukawa type theory is equivalent to a purely 
fcrmionic model. For this purpose we perform the Gaus- 
sian functional integration over the <fi field. Expressed 
only in terms of fermions our model contains now a four- 
fermion interaction with a momentum dependent cou- 
pling A(Qi, Q2, Q3, Qi). The interaction has the explicit 
form (Q 4 = Qi + Q 2 - Q 3 ) [H 



c I 



[ X 4> 



hi 



27ri(m - n 4 )T + (q\ - q^f/AM + D A - 2a 



}■ (3) 



In the pointlike limit the momentum dependence can be 
neglected and eq. J3J is replaced by the "local interaction 
approximation" 



A = A„ 



va — 2(7 



(4) 



Later we will argue that the pointlike limit can be estab- 
lished by h,/, — > 00, D oc (with D the UV renormalized 
version of Da, cf. 0] and sect. IIII A|l . This settles the 
connection between broad Feshbach resonances and the 
pointlike limit, for which a purely fermionic description 
as discussed by Strinati et al. [l3. 120L I2H I22I l23l | is ap- 
propriate. 

The action has a global E7(l) symmetry. This 

implies the conservation of the particle number N — 
J d 3 x n, where the total number density of atoms n in- 
cludes the unbound atoms and those bound in molecules 
or forming condensates. It obeys 



n F + n B = (^ip) + 2(cf>*4>) 
n F + 2fi M + n C - 



(5) 



In the second line, we have expressed the disconnected 
correlation function (4>*4>) by the connected Green's func- 
tion (4>*4>) c and the expectation value (0). These terms 



can be interpreted |18| as density contributions from open 
channel atoms {rip = (ip'ip)) and closed channel atoms in 
uncondensed molecules (um = {<t>*<t>) c) as well as a con- 
tribution from closed channel atoms in the condensate 
(n c = 2(0*) (<£)). We emphasize that cq. © is no "ad 
hoc" assumption - it directly follows from the microscopic 
formulation of the functional integral. 

For a homogeneous setting the condensate is constant, 
(4>(x)) = 4>q. It can be related 18] to the expectation 
value of a composite di-atom operator 



2va - 4cr 



(6) 



This demonstrates directly that our formulation in terms 
of a path integral makes no difference between a "conden- 
sate of molecules" </> and a "condensate of atom pairs" 
(ijj T eil>} - they are simply related by a multiplicative con- 
stant. Hence, the universality of the condensation phe- 
nomenon itself - the fact that we do not need to specify 
a precise microscopic mechanism as e.g. the formation of 
Cooper pairs - is particularly visible. 

The field equations which determine 4>q and n can be 
computed from the effective action T[a, <fi], which depends 
on a and the expectation value cf> = (</>). As usual, it 
obtains by a Legendre transform, eliminating the source 
j(j,(x) in favor of <j>(x), 



T[<f>] = - lnZufoM] + / dx{jl<t> + U4>*)- (7) 

For homogeneous a, </>, it is appropriate to consider the 
effective potential U = (T/V)T. To one loop order it is 
composed of a classical piece and a fluctuation contribu- 
tion, 



U = U {cl) +Ux. 



(8) 



In our case, the classical contribution reads 

= (D A - 2a)^*l (9) 

The one loop contribution is composed of a piece from 
fermion (U\) and boson (u[ B ^) fluctuations. The ex- 
plicit expression for the one loop part depends on the 
detailed approximation. It will be discussed in sect. IVII 
while we will first display some general features and re- 
sults of our approach. 

The homogeneous field equations read 



dU 



dU 



n, 



= 0. 



(10) 



The condensation of atoms pairs or molecules is signalled 
by a non- vanishing expectation value ((f>) = 0o, given 
by the location of the minimum of U. The associated 
spontaneous breaking of the global continuous symmetry 
of phase rotations of ip and (f> (related to the conservation 
of the number of atoms) induces a massless Goldstone 
boson. This is the origin of superfluidity. 
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III. UNIVERSALITY 



Additive Renormalization 



The microscopic action depends explicitly on three 
parameters V\, h<f,, X^. A fourth parameter is introduced 
implicitly by the ultraviolet cutoff A for the momentum 
integration in the fluctuation effects. (Besides this, the 
results will depend on the thermodynamic variables T 
and a - see [lg for the relation between the effective 
chemical potential a and the chemical potential //.) The 
renormalization procedure eliminates the cutoff depen- 
dence and trades these "bare" parameters for renormal- 
ized ones which can subsequently be related to observ- 
ables of the concrete atomic system (cf. sect. |X|) . 

It turns out that it is sufficient to perform the ultravi- 
olet renormalization in two instances: The first one con- 
cerns the connected correlation functions (density contri- 
butions) (tp^tp) an d (^*4>)c- Carefully relating these field 
expectation values to their counterparts in the operator 
formalism shows that the div erg ences cx A 3 are due to a 
zero point shift of the density i8j . In our setting the con- 
tributions of the fermionic fluctuations are cancelled by 
the bosonic (molecule) fluctuations and no renormaliza- 
tion is necessary in this respect. The second divergence 
for A — > oo is associated to the bare detuning T>a- It can 
be absorbed by using the renormalized detuning 



2tt 2 ■ 



This prescription corresponds to a renormalization of the 
effective atom interaction strength A Q} |18| . 

Expressed in terms of v and h$ the effective potential 
becomes very insensitive to the microscopic physics, i.e. 
the value of the cutoff A. Without much loss of accuracy 
we can take the limit A — > oo for the computation of U. 



B. Dimensionless Parameters 

The action still involves the atom mass M which 
is much larger than all other characteristic mass scales. 
By the choice of appropriate units the nonrelativistic 
physics becomes independent of M. We therefore in- 
troduce a "Fermi momentum" kp and a "Fermi energy" 
ep = k F /2M. The inverse of the Fermi momentum 
should be associated with the most important charac- 
teristic length scale in our problem. Roughly speaking, 
it corresponds to the average distance d between two 
unbound atoms or molecules. Nevertheless, the precise 
value of kp is not relevant - it just sets the units and we 
could take a fiducial value fcp = leV. For a homogenous 
situation it is convenient, however, to define hp by the 
density, 



Next we rescale coordinates and fields according to 

x = kpx, f = epr, T~T/ep, q = q/kp, (13) 
tp = k~ 3/2 tp, <p = k~ 3/2 (f>, a = a/e F . 

This yields the scaling form of the action 



Sr = dx 



^(B T -A-a)i> 



-4>*{B T 



1a 

2 



D A - 2a) 



(14) 



Accordingly, all quantities derived from the partition 
function can be brought to a scaling form and can only 
depend on the dimensionless parameters {Da = v + 

Kiy^kp)) 



v/e F , 



2Mk 



-1/2 



V \^ = 2Mk F X 4> . (15) 



This scaling form shows a first important aspect of uni- 
versality. All computations can be performed at a fidu- 
cial kp — leV. Contact to a concrete physical situation is 
only made at the end by rescaling the results by using the 
value of kp appropriate to the density of a given experi- 
ment. In our approximation for a homogeneous situation 
all dimensionless physical observables can be expressed 
(11) in terms of the dimensionless variables v, h^,, X^ and T\ 
The scale is then introduced by n. 

The parameters in the action (|14fl are microscopic pa- 
rameters VA-, /i0,Aj \l>,K- Similar to the renormalization 
of v also the couplings and A^ will be affected by 
renormalization effects. For a comparison with experi- 
ment we have to relate the "bare couplings" h^Ai A^a 
to renormalized couplings, as for example the coupling 
h^fi which appears in the molecular binding energy in 
vacuum of the scattering of two isolated atoms 3 . We 
will argue below that in terms of the macroscopic rele- 
vant parameters only two parameters remain which can 
be associated to the renormalized couplings v and h^. 
The macroscopic value of A,/, becomes computable as a 
partial fixed point in terms of v and hj,, at least suffi- 
ciently close to the resonance and for sufficiently large h$. 
Furthermore, for the broad resonance limit of a pointlike 
interaction we may absorb X^ by a Fierz transformation 
into a redefinition of h \ such that X^ Q remains invari- 
ant, i.e. h% — ► TA — X^va (with \2a\ <C Pa)- In sect. [X] 

we will relate v and h<f, to the parameters of experiments 
with gases of 6 Li and 40 K. 



k F /(3n 2 ) 



(12) 



3 In order not to overload the notation we often will not specify 
directly the indices like /i^ a or ^</>,o ■ Where it becomes necessary 
we describe the status of more precisely 
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C. Concentration 

The space of model parameters and thermodynamic 
variables for the crossover problem is spanned by v,f& 

and T. However, it turns out that another set of param- 
eters, c, and T, describes the system more efficiently, 
revealing further aspects of universality more clearly. Us- 
ing the relation between a four-fcrmion coupling A and a 
scattering length a , 



A 



Ana 



(16) 



we define an effective in-medium resonant scattering 
length a or, in the scaling form, 



1 

akp 



8tt(z>- 2a) 



(17) 



We will refer to c as "concentration". Indeed, with the 
inverse Fermi momentum setting the scale for the av- 
erage interparticle spacing d, the concentration c is a 
measure for the ratio between scattering length a and 
average distance, c oc a/d. It is therefore of no surprise 
that in the case c^ 1 < 1, where the scattering length 
exceeds the typical interatomic distance, many-body ef- 
fects play an essential role. Note that through a, this in- 
medium scattering length depends on the density. It will, 
in general, differ from the "reduced scattering length" 
<ir = a(<7 = 0), which is defined by 



1 



ankf 



IT' 



(18) 
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FIG. 1: Crossover phase diagram. The dependence of the 
critical temperature T c = T c /eF on the inverse concentra- 
tion c" 1 = (afcF)" 1 is shown for the broad resonance limit 

h 



4> 



oo (solid and short-dashed line) and for the narrow : 



onance limit — > (dashed-dotted line). We also indicate 
the standard BEC (dashed horizontal line) and BCS (dashed 
rising line) limits which do not depend on the choice of the 
Yukawa coupling. For the broad resonance limit we plot two 
different approximations. 



In summary we have now an effective low energy for- 
mulation where neither M nor A enter anymore. Every- 
thing is expressed in terms of kp and three dimensionlcss 
parameters, namely c, T and h^. 



IV. ENHANCED UNIVERSALITY 



For broad Feshbach resonances the concentration is re- 
lated to the full scattering length a(B) by c — a{B)kp (cf. 
sect. |Xj. Here the dependence of a(B) = a res (B) + a,b g 
on the magnetic field B involves both a resonant piece 
a res (B) that diverges at the Feshbach resonance B — > Bo, 
and an approximately B - independent background scat- 
tering length cibg. 

The impact of the choice of variables c~ 1 ,h^ ) ,T is 

twofold. First, it absorbs a good deal of the - de- 
pendence in the final results. In the language of crit- 
ical phenomena, the dimensionless temperature T and 
the inverse concentration c _1 are relevant parameters, 
whereas /i^ is marginal. A more systematic discussion of 
the role of will be given in the next section. Second, 
the choice of c or a relates our formulation based on par- 
tial bosonization more closely to an interacting fermion 
system, where the effective scattering length a is obvi- 
ously a crucial parameter. 

For a small concentration |c| the gas is dilute in the 
sense that scattering can be treated as a perturbation. In 
this range mean field theory is expected to work reason- 
ably well. On the other hand, for large |c| the scattering 
length exceeds the average distance between two atoms 
and collective effects play a crucial role. 



Enhanced universality occurs in situations where one 
of the two parameters c, h<j, becomes irrelevant (or fixed). 
This concerns the limits of narrow and broad Fesh- 
bach resonances ( "narrow resonance limit" , — > and 
"broad resonance limit", — > oo), the BCS and BEC 
regimes, c — > 0, and the scaling limit, c — » oo. For a 
given T the crossover for narrow and broad Feshbach 
resonances can be described by a single parameter c _1 . 
Results for the BCS, BEC and scaling limit can only de- 
pend on hrf, and become parameter free for narrow or 
broad resonances. 

We can qualitatively discuss the limits of enhanced uni- 
versality taking the phase diagram in fig. ^as an exam- 
ple. The values of c _1 and determine the phase di- 
agram completely, i.e. the critical temperature depends 
only on these two parameters, T c — T c (c _1 , h^). We dis- 
play the critical temperature T c as a function of c _1 and 
compare the limits — > (dashed) and — > oo (two 
approximations, solid and short dashed). (The underly- 
ing approximation schemes are discussed in sect. IV1I ) 
For intermediate values of h^, a monotonic change be- 
tween the two limits is found, cf. also fig. [21 Though the 
impact of the marginal parameter is moderate, the 
precise value of hs influences the details of the crossover 
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and will be important for precision estimates for concrete 
physical systems, in particular when narrower resonances 
than those presently investigated in 6 Li or 40 K will be ex- 
plored. 

Fig. ^dearly reveals an h$ - independent approach to 
the limiting BCS and BEC regimes: for | c 1 1 — > oo, h$ 
becomes an irrelevant parameter and the system can be 
described in terms of the concentration only. For com- 
parison, we have also plotted the extrapolated standard 
BCS (dashed, rising line) and BEC (dashed, horizontal 
line) results. 

The narrow resonance limit corresponds to an exact so- 
lution of the many-body problem. It requires — » and 
additionally a vanishing four-fermion background inter- 
action, but is free of further approximations. The broad 
resonance limit, instead, corresponds to a strongly in- 
teracting field theory, whose approximate solution still 
involves quantitative uncertainties, in particular close to 
the critical temperature (cf. sect. IVI|) . This is reflected 
by the two critical lines obtained in different approxima- 
tions in the broad resonance limit. The difference be- 
tween the solid line and the short dashed line gives an 
idea of the remaining uncertainty in the treatment of the 
molecule fluctuations. For the solid line our result at the 
resonance T c = 0.292 is in agreement with the theoret- 
ical value obtained in |4^, |4g and compatible with the 
measurements reported in ^ij T c = 0.31 ± 0.04. Includ- 
ing the molecule fluctuations in the gap equation (short 
dashed line) yields T c (c _1 = 0) = 0.255. The maximum 
of T c (c _1 ) has a similar value for both approximations. 

We find a new form of crossover from the narrow to 
broad resonance limit in dependence on the Yukawa or 
Feshbach coupling h^. This is shown in fig. |2 where 
we plot T c (fig. [21 (a)) as a function of h$ for c _1 = 0. 
(The curves for other values of large |c| are similar.) One 
clearly sees a smooth interpolation between the regimes. 
Within the broad resonance regime the precise value of 
/i0 is not relevant for T c or, more generally, for all quanti- 
ties not involving explicitly the microscopic constituents 
("bare fields"). Renormalization effects for h& even as 
large as a factor of ten do not matter. In contrast, in the 
"crossover regime" of intermediate h$ ~ 10 the physical 
results depend on the value of h^. Now renormalization 
effects must be taken into account carefully for a reli- 
able computation. Furthermore, we will learn in sect. 
El that observations related to "microscopic observables" 
like the fraction of closed channel atoms (or microscopic 
molecules) can depend strongly on even in the broad 
resonance regime. We also show in fig. [21 (b) the gradient 
coefficient for T = T c , cT l = - this quantity will be 
defined in the next section. 

More details on the different limits of enhanced uni- 
versality will be given in sect. IVIIII after the discussion 
of the necessary formalism in the next sections. 
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2: Enhanced universality for large and small h$. For 
c " = we plot the dependence on of (a) the dimension- 
less critical temperature T c and (b) the gradient coefficient 
A$ for T = T c (cf. sect. [VJ. For small h</, < 1, a sta- 
ble universal narrow resonance limit is approached. For large 
h$ we find a very pronounced insensitivity of T c and A^ to 
the precise value of h$ - note that we plot on a logarithmic 
scale. The "crossover" regime interpolates smoothly between 
the universal limits. 



V. DRESSED AND BARE MOLECULES: THE 
WAVE FUNCTION RENORMALIZATION 2$ 

For an understanding of the crossover problem, it 
is crucial to distinguish bare and dressed molecules 
ItH |48| . In our interpretation, the bare or "closed chan- 
nel" molecules are associated to a true nonlocality in the 
microscopic interaction |[3J. Their overall role for the 
thermodynamics should be subdominant for broad reso- 
nances or almost pointlike microscopic interactions. The 
dressed molecules can instead be effective degrees of free- 
dom, generated dynamically by fluctuation effects. They 
correspond to a renormalized bosonic field. 



A. Derivative Expansion 

For a homogeneous situations or weakly varying trap 
potentials a derivative expansion of the effective action 
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becomes appropriate (with (f) 
of sources j^) 



(4>) = k F 3 ^ 2 (j) in presence 



dx{u(4>) + AtpV^Vj) + Z^*d T 4> + ...}. (19) 



By analytic continuation from euclidean time t to 
Minkowski time ("real time") we can derive a time 
evolution equation for <fi from the variational principle 
ST/S(j> = 0. The coefficient of the term with the time 
derivative, Z^, plays the role of a wave function renor- 
malization. We therefore introduce renormalized fields 
for bosonic quasiparticles or "dressed molecules" and ap- 
propriately renormalized couplings 



= z l /2 & p = 4>*4>, M = 



v=^T (20) 



such that 



r[(f>} = J dZ{u{cj))+A^^V (j) + ^d T (j) +...}. (21) 

The effective potential can be decomposed into a clas- 
sical contribution contribution from the fermion 

fluctuations ii\ and a contribution from the molecule 



fluctuations 



u = 



(22) 



[v - 2a)4>*4) = {y- 2a/Z <t> )p. 



The a - dependence of the different pieces is closely re- 
lated to the different contributions to the total density 
of atoms [l8| . Indeed, the pressure of the system obeys 
p(T) = -\Uo((T,T) - Uo(a = T = 0)), where U (a) de- 
notes the value of the effective potential at its minimum 
(with respect to <fi) . A thermodynamic relation associates 
the particle density with the derivative of p with respect 
to the chemical potential. 



B. Contributions to the atom density 

The wave function renormalization Z$ plays a crucial 
role for the precise understanding of the different con- 
tributions to the total atom density, from condensed or 
uncondensed atoms, molecules or unbound atoms, open 
channel and closed channel atoms. For example, the 
number density of uncondensed dressed molecules um 
can be computed from the propagator for the renormal- 
ized field <p. It is related to the a - derivative of [Tsl ] 
evaluated at fixed v (not fixed c _1 !) 
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FIG. 3: Fraction of closed channel molecules = Qc +f!jvf, 
compared to experimental values for T = and fci? = 
0.493eV=250nK. On the BEC side, we find that our results 
are quite insensitive to the precise choice of &f- The strongly 
interacting region c _1 < 1 is indicated by vertical lines, where 
the center line denotes the position of the resonance. For 
B > Bo the solid and long dashed-dotted lines reflect the 
uncertainty from the "Fierz ambiguity" discussed in app. \B\ 
For B < Bq the dashed line omits the renormalization effects 
discussed in sect. DO and app. 



molecules. Similarly, the number of bare or microscopic 
molecules um corresponds to the propagator for <f> such 
that 



n M = Zd,nM, ^ 



M 



(24) 



Now VLm counts the fraction of closed channel atoms (not 
in the condensate) as defined by the microscopic molec- 
ular Feshbach state. For large Z^, the number of dressed 
molecules um can be substantial even if the number of 
microscopic molecules um is negligibly small. Of course, 
in this case the dressed molecules are dominantly com- 
posed of open channel atoms and not of closed channel 
molecules. 

Similarly, the condensate fraction Qc = nc/ n mea- 
sures the fraction of atoms in the condensate 



3^2 = fcT = 2 ^o = 2 Po = -Z, 



(25) 



It is related to nc {S} by nc = Z^hc- Thus the fraction 
of closed channel atoms in the condensate obeys f2c = 
n c/ n — ^c/^4> an d is tiny for large Za,. The remaining 
number of (dressed) unbound atoms obeys np — n—riM— 
nc- Correspondingly, Qp = np jn denotes the fraction of 
atoms that are neither condensed nor bound in dressed 
molecules. 



hp 



,-.(») 



ldu\ 
2 



da 



67?" 



(23) 



Here </>o corresponds to the minimum of u and Qm = 
2um I n counts the fraction of atoms bound in the dressed 



C. Computation of Z$ 

The wave function renormalization Z$ can be com- 
puted from the dependence of the molecule propagator 
on the (real) frequency |18|. We include in this paper 
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only the contribution from fermion fluctuations. In this 
section we employ the relation [l^ 



1 d 3 u[ 



rS F ) 



2 dad<t>*d(t> 



(26) 



which relates to the dependence of the inverse 
molecule propagator on the effective chemical potential 
a . (In the app. [Owe discuss a definition of Z^ based on 
the explicit frequency dependence of the inverse propa- 
gator.) In eq. (12(:i[) is computed from a fermion loop 
in the background of constant 4> and a. It reads 



~ (F) i ~ 



= -2T 



d 3 q 

(27T) 



In (e^ 7 + e~ 7 *~ 7 ) 



ATq 2 



(27) 



The limit T — > 0, a — > is infrared divergent. This di- 
vergence is cured by a nonzero momentum q in the prop- 
agator: For T = a = the correction oc h^q 2 )^ 1 / 2 
diverges only for q 2 — > 0. In the vacuum the physics of 
the "dressed molecules" has therefore to be handled with 
care. We will see in sect. IIXI that the case a — ► ap- 
plies to a situation where the molecules correspond to an 
unstable resonant state with positive binding energy. In 
contrast, stable molecules correspond to a nonzero neg- 
ative value of <7 without infrared divergence in Z^. Fur- 
thermore, for nonzero density a condensate occurs for 
T = and Z$ is regularized by the nonvanishing value 
of f in eq. 



D. Renormalized Yukawa coupling 



with 



70 
7 



2 s 1/2 
+ f 



(28) 



70(0 = 0) 

v = h%$% 

The relation l|26Jl yields 
hi r d 3 q 



2T X 



Z A = 1 



16T 2 J (2tt) 



770 3 [ tanh 70 - 70 cosh 2 70 ] 

(29) 



and agrees indeed with the definition from the frequency 
dependence of the propagator in app. |0 

In the symmetric phase Z$ can be expressed in terms 
of the dimensionful parameters as 

z = i + lM 3 / 2 ^r%( s ), 



y 



{ tanh(z/ — s) g — ^ 7 } > 



cosh (y 2 — s) 



2T 2T 
where -ff(s) has the limit 

i?(s -> -00) 



|(-)- 1/2 - 



(30) 



(31) 



4 The appropriate choice of the Yukawa coupling relating r and 
is a subtle issue. We follow here the concept of a systematic 



renormalization group improved perturbation theory 49] where 
/i0 is a a and T dependent renormalized coupling which is evalu- 
ated at the momentum scale which dominates the integral 1271 . 
This prescription differs from the Schwinger-Dyson equations for 
the fermion contributions to the <f> vertices which follow from the 
- derivatives of eq. 1271 . In the SDE - approach one factor of 
would correspond to a microscopic or bare coupling fi.^ ; A- 



We observe that for large h$ the wave function renor- 
malization diverges Z$ oc hi. We will find that the total 
fraction of dressed molecules ils — Hm + is a quan- 
tity of the order one in the crossover region. We con- 
clude that in the broad resonance limit the total number 



of bare molecules = i! 



9. ( 



Qb/Z^ 



becomes negligibly small. As argued in sect. IIVI the 
microscopic molecules can be omitted for the thermody- 
namics and the broad resonance limit is well described 
by a pointlike interaction for the fermionic atoms. Nev- 
ertheless, the renormalized bosonic field <f> remains a very 
useful concept. We will see that in the limit h<f, — > 00 the 
appropriately renormalized couplings remain finite and 
become therefore independent of the precise value of h^. 
This holds, in particular, for the renormalized Yukawa 
coupling 5 



z- 1/2 h 



and the gap 

r = h 2 J*4> = (hl/Z^Wt = hi, 



(32) 



(33) 



The fact that these quantities are effectively - inde- 
pendent has a simple but important implication regard- 
ing the renormalization effects for hj, itself: Even if the 
actual value of is changed by, say, an order of magni- 
tude, this will not affect these parameters. All quantities 
that can be expressed in terms of the renormalized cou- 
plings become very insensitive to the value of h^, as seen 
in fig. 



5 The renormalization of the Yukawa coupling is done in two steps. 
A first renormalization step computes a renormalized coupling 
h^, as extracted from binding energy and scattering length in 
vacuum, e.g. hc/yfi in sect. |X| A second step 1321 takes the wave 
function renormalization into account. 
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E. Bosonic quasiparticles 

We have computed the gradient coefficient in [Isj . 
In fig. El(b) we have displayed the renormalized gradient 
coefficient Aa, at the resonance, c^ 1 — 0. It switches from 
0.5 for small to a value close to one for large h$. On 
the other hand, deep in the BEC limit we find A^ = 0.5 
both for small and large h$. Indeed, the BEC limit is 
well described by pointlike bosonic quasiparticles, even if 
the microscopic molecules become irrelevant for h<f, — > oo. 
The limit 



lira lira A^ 



(34) 



states that the dressed molecules always have mass 2M. 
Furthermore, in the BEC limit the number density of 
atoms (<j)R = Z 



1/2, 



n = 2n M + n c = 2({fa$ R ) e + (4>* r )(4>r)) (35) 

turns out to be twice the number of dressed molecules 
(including the ones in the condensate) according to the 
standard Bogoliubov formula. (For more details, we refer 
to ^| . A similar result has been obtained by Strinati el 
al. in a purely fermionic setting |2£|[2^.) In the BEC 
limit all quantities that can be expressed in terms of the 
renormalized couplings can equally well be described by 
a simple model involving only bosonic effective degrees of 
freedom. This holds both for the case where microscopic 
molecules dominate and for the case of a purely point- 
like interaction for fermionic atoms without microscopic 
molecules. In this case an effective bosonic theory looses 
memory of its microscopic origin. 



F. Fraction of microscopic molecules 

Our calculations relate the "macroscopic observables" 
and the renormalized parameters to the microphysical 
interactions. In view of the large universality for the 
experimentally tested broad resonances one may won- 
der if these calculations can be checked experimentally 
or if only universal relations will remain at the end. A 
first possible check of the detailed relation between the 
microphysics and macrophysics concerns the relation of 
the universal concentration parameter c with the detailed 
experimental setting, in particular its dependence on the 
magnetic field B. This will be discussed in sect. [X] A 
second access is provided by tests that directly observe 
"microscopic quantities" as, for example, the fraction of 
microscopic molecules. 

Indeed, even for a broad Feshbach resonance the mi- 
croscopic molecules exist and one can measure their num- 
ber. The laser probe in the experiment by Partridge et 
al. triggers resonant transitions from closed channel 
microscopic molecules to open channel atoms. It there- 
fore directly couples to the total number of microscopic 



molecules ((j>*(j>}. The quantity Z in [j| can be identified 
with the fraction of bare molecules 



Ci B = n M + = («m + ^c)/Z^. 



(36) 



We have computed this quantity as a function of c _1 and 
relate c _1 to the magnetic field in sect. [X] The result 
is shown in fig. The comparison with the experimen- 
tal values is very good. Since Q-m + is of the order 
one the dominant effect actually checks our computations 
of Z$, Obviously, Ob depends explicitly on the Yukawa 
coupling, £lg oc /z.0 2 , in distinction to the universal quan- 
tities for a broad Feshbach resonance. For T = and 
n = the Yukawa coupling h^^ can be extracted from 
properties of single atoms and molecules (cf. sect. |Xj 
and one finds for 6 Li hfj} = 610 (k F = leV). However, 
even for T — (as appropriate for the experiment [j|) the 
renormalization effects induce an effective dependence of 
htf, on c _1 as shown in app. El Omitting these effects 
gives substantial disagreement with the measured values, 
as shown by the dashed curve in fig. [3] for B < Bq. We 
consider fig. [21 as an important quantitative check of our 
method 6 that we will describe in more detail below. 



VI. EFFECTIVE POTENTIAL 

In our approximation all relevant information is con- 
tained in the various contributions to the effective po- 
tential u^ cl \u^\u^ and in A^Z$. These quantities 
depend on T, a, c _1 and h$ where a can be fixed by the 
first equation (fTTTf) 



With f 



h\p, (8ttc)- 



du _ 1 
da 37r 2 

1 = -hr 2 (v - 2a) and 



8nc 



dp,* ^ 



\f,a) 



-(B) 



(P) 



(37) 



(38) 



we see that for fixed c the explicit dependence of £t(f) 
on hfj, enters only through the contribution from the 
molecule fluctuations u\ . The mean field theory (MFT) 

( B) 

only includes the fermionic fluctuations and omits u\ . 
According to eq. I|23|) standard MFT also omits the den- 
sity of unbound dressed molecules, i.e. Om = 0. With 



Of,o — 

n c = 



-3tt 



-37T* 



da ' 



6ir 2 po 67r 2 fo 



da 



(39) 



6 Another quantitative check concerns the agreement or o ur resul ts 
for T = llSl with Quantum Monte Carlo simulations l5Ctl5ll . 



11 



eq. (|37|l is equivalent to Of.o + = 1- For small (or 
zero) Clc the dependence on hs becomes negligible. The 
MFT equations simply reduce to 



3^ 



, du 



r,(1 



da 



= 1- 0, 



- (40) 



with ; given by eq. 1)27(1. 

A first step beyond MFT includes the contribution of 
in eq. (|37|1 and omits it in the field equation du/d(j)- 
This changes the first equation (|40[1 to 




Qf o = — 37r 



,du\ 



= 1 - o M - o 



(41) 



The presence of a substantial Oa/ on the BEC side re- 
duces n^o and therefore changes the value of a. This 
effect is already sufficient to produce a crossover from 
the BCS to the BEC regime. We therefore need an esti- 
mate of Ojyf. 

According to eq. (J5J the bare molecule density is re- 
lated to the exact propagator of the bosonic field 



= 9 tr E 



d 3 q 
(2tt 3 ) 



G,p(q,n) - h B . 



(42) 



with ns an additive renormalization .18]. Here we work 
in a basis with dimcnsionlcss renormalized real fields 
01,02 with_0 = (0i+i0 2 )/\/2 andp = 0*0 = (0?+0|)/2. 
Therefore = "P^ 1 is a 2 x 2 matrix over which tr takes 
the trace. In terms of renormalized dimcnsionless quan- 
tities, V& = P^/iZ^ep), this yields for the density of 
dressed molecules 



hp 



T 



-tr 



E 



(2tt 3 )' 



p 

F 



(43) 



This relation is exact. It involves the exact inverse propa- 
gator P^, which is directly related to the second functional 
derivative of the effective action (|21|l with respect to 0, 



S 2 T 



54>*a{z)Hb{x'y 



(44) 



We will employ a Taylor expansion of the effective po- 
tential u around its minimum 



rn 4>P+ \ X <PP 2 + 



SYM 
SSB 



(45) 



where SYM and SSB denote the symmetric and super- 
fluid (spontaneously symmetry- broken) phases. Without 
loss of generality we also choose 0o real such that 0i cor- 
responds to the "radial mode" and 02 to the "Goldstone 
mode" . In the approximation (|21|l this yields 01 



A^q 



r,2 



m 4 



2A rf 



OJn 
i ™2 



, A A q 2 + m 



(46) 
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FIG. 4: Temperature dependence of the bosonic mass term 
rr& for c _1 =0,h$ — > oo. The vanishing of the mass term 
is approached continuously (solid line) . We also plot the sepa- 
rate contributions from "mean field" (classical plus fermions, 
dashed) and boson fluctuations (dashed-dotted). For T > T c 
the latter become unimportant away from the phase transi- 
tion. 



and we note the appearance of the Matsubara frequencies 
in the off diagonal entries. Here we use a formulation 
that can be used in both the normal (SYM) and the 
superfluid (SSB) phase. We will specialize to these cases 
in the following subsections. The mass matrix 



d 2 u 



2X 4 




m 



(47) 



involves both the fermionic and bosonic fluctuation cor- 
rections to u. At the phase transition (jrA) a b vanishes 
and we obtain 



47T^ 



(48) 



The field equation for the renormalized field expecta- 
tion value can be expressed as 



= 0, 



In the superfluid phase 



iii- 



m l = 



du 
dp 



(49) 



7^ 0, po 7^ 0) one infers 
0. Therefore the mass matrix <|47|) has a zero eigen- 
value and 02 can be identified with the massless Gold- 
stone boson. The presence of a "massless mode" with 
infinite correlation length is the central ingredient for su- 
perfluidity. For the symmetric phase the two eigenvalues 
of (mV) a b are equal since po = 0. 

In principle, the inverse propagator H4fi|l involves a mo- 
mentum dependent vertex X^(q). Including the molecule 
fluctuations \$(q) vanishes for q = 0. As advocated 
in |l8j we take this difficulty into account by using in 
the density equation 143|) only the fermionic contribu- 



tion Xa 



A 



(F) 



d 2 u[ F) /dp 2 \ po . To reduce the error 
connected to this approximation one may either com- 
pute the momentum dependence of the four-boson ver- 
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(T - T c )/T c 

FIG. 5: Molecular correlation length £m as a function of the 
reduced temperature (T — T c )/T c . We consider three regimes: 



crossover (c = 0, solid), BCS (c 



-1.5, dashed), BEC 



(c 1 = 1.5, dashed-dotted) and the broad resonance limit 



tex or, perhaps more promi sing , use functional renormal- 
ization group methods EH ■ The approximation 

(F) 

= is, of course, consistent with our "intermediate 
approximation" where the field equation for (f> neglects 
the bosonic fluctuation contributions. 

A more consistent inclusion of the molecule fluctua- 
tions also takes them into account in the field equation 
for 4>, modifying also the second equation i|40|) . We need 
to evaluate the bosonic mass term 



(B) 



dp 



(F)2 



du[ B) 

dp ■ 



(50) 



While Ui is given by eq. I|27|l we still need to determine 

du\ I dp. For this purpose we compute the lowest or- 
der "gap" or Schwinger-Dyson equation for the molecule 
propagator. In the context of a perturbation expansion 
the gap equation amounts to a self-consistent resumma- 
tion of a series of diagrams involving bosonic fluctuations 
beyond the Gaussian (one loop) level. The explicit form 
of the gap equations for the normal and superfluid phases 
as well as the phase boundary is displayed in app. [Bj De- 
tails of their derivation and limitations of validity can be 
found in [l8|. 



VII. CORRELATION LENGTH AND 
SCATTERING LENGTH FOR MOLECULES 

In the symmetric phase the static effective molecule 
propagator reads (for zero Matsubara frequency) 



1 



Aj>q 2 



(51) 




BEC" 



- 1 - 0.5 0.5 1 1.5 



FIG. 6: Effective potential in the crossover region, T = 
0.33, htj, — * oo . The solid line expresses the quartic coupling 
A</> in terms of the renormalized molecule scattering length ob- 
tained with clmIcf = A^/(47r). The short dashed line neglects 
the molecule fluctuations and shows ar^pkF = A^ F ' / '(47r). 
Large renormalization effects are observed in the crossover 
regime. In the BEC limit, the bosons are only weakly inter- 
acting. In addition, we plot the renormalized mass term for 
the molecule field rri^ (dashed-dotted) and the renormalized 
gradient coefficient A$ (long dashed). 



We approximate here A$ as a constant. The Fourier 
transform of G$ decays for large distances r exponentially 



G^(r) ~ exp( 
with correlation length 



<mt) 



M 



V A 



1/2 



(52) 



(53) 



(More precisely, £Tf corresponds to the location of the 
pole of the propagator for negative q 2 and therefore 
should be evaluated at the location of the pole instead of 
q 2 = as approximated here.) At the critical tempera- 
ture T c the correlation length £ diverges since to^ = 0. 

In the superfluid phase the propagators for the radial 
and Goldstone modes differ. The correlation length for 
the Goldstone fluctuations (<fo) is infinite according to 
eq. (|49|) . In the "radial direction" (<fii) the effective 
mass term for the fluctuations around the minimum reads 
2A0po. H vanishes (at zero momentum) since our approx- 
imation based on Schwinger-Dyson equations [18| yields 
A0 = for the deep infrared value of the molecule cou- 
pling. Hence we obtain in our approximation a diverging 
correlation length for the radial mode. Nevertheless, the 
pole is actually not located at q 2 = - using a momentum 
dependent A^(g 2 ) in m 2 ^ or taking this effect into account 
in the form of A^q 2 ) would be a better approximation. 
In this paper we will not discuss further the effects that 
render the radial correlation length finite for T <T C . 

The dimensionless molecular correlation length in the 
symmetric phase 



(54) 
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is expressed in terms of the renormalized quantities 
and mi 



= A 



A? 



A 



3tt 2 



n 



M 



}. (55) 



For the computation of A^ we refer to 01 ■ We plot the 
renormalized boson mass, displayed in curly brackets in 
eq. 1)55(1 . in fig. Additionally, we indicate the separate 
contributions from fermion (dashed) and boson fluctua- 
tions (dashed-dotted) . Boson fluctuations are important 
close to the critical temperature. 

The condition £ M 2 = defines the critical tempera- 
ture where the molecule correlation length diverges. Fig. 
displays £m in the BEC, crossover and BCS regime. 
We observe that the molecular correlation length re- 
mains larger than the average distance between atoms 
or molecules ~ fc^ 1 even rather far away from the crit- 
ical temperature. This feature is quite independent of 
the BEC, BCS and crossover regime. It indicates the 
importance of collective phenomena. 

For £m 3> 1 the correlation length exceeds by far the 
average distance between two atoms. We do not expect 
the loop expansion to remain accurate in this limit. This 
is the typical range for the universal critical behavior 
near a second order phase transition which is known to 
be poorly described by MFT or a loop expansion. A 
proper renormalization group framework is needed for 
£m S> 1- In contrast, for £m <C 1 the molecules are 
essentially uncorrelated. 

Let us define an effective molecule scattering length 
a m for the dressed molecules in analogy to eq. I|16|) (with 
"classical" boson mass 2M for composite particles) 



Ana 



M 



2M 



such that 



A,* = 



a Mk F = a$kF = -*^-. 

A second possible criterion for the validity of the loop ex- 
pansion is the smallness of the molecular self interaction 
A^. Typically, a one loop expression becomes question- 
able for ajvffcj? 1. 

We plot clm^f and a^kp in fig. E|as a function of cr 1 
and in fig. as a function of temperature. Both plots 

( F) 

are for the broad resonance limit. While a M kp grows 
large near T c the "full" molecule scattering length aukF 
goes to zero. The reader should be warned, however, 
that for T — > T c the momentum dependence of \<f,(q) 
becomes crucial - for T = T c one has A^(q) cx \fq* such 
that A^ vanishes only for q — 0. At T c the interaction 
between the dressed molecules cannot be approximated 
by a pointlike interaction. 

Finally, we comment on the ratio of fermionic and 
molecular scattering length in the two-body limit. The 



2Mk ? 



(56) 




(T - T c )/T c 

FIG. 7: Scattering length for molecules aMkF as a function of 
the reduced temperature (T — T c )/T c in the crossover regime 
(c _1 = 0, htf, — > oo,T c — 0.255). Boson fluctuations renormal- 
ize the bosonic scattering length down to zero in the vicinity 
of the critical temperature. This may be compared to the 
increase of the fermion fluctuation induced part a/y kp. 



diag rammatic approach of |l9j | yields ajv//a ~ 0.75. In 
|55| , the ratio was extracted from an analysis of the quan- 
tum mechanical four-body problem in the BEC limit, es- 
tablishing the result am /a w 0.6. This was reproduced 
in a Quantum Monte Carlo simulation in |5lj . 

We discuss a procedure to project on the low den- 
sity limit in sect. IIXI For the fermionic contribution to 
the molecular scattering length, this yields in the BEC 
regime 



(57) 



This is the Born approximation. In the low density limit, 
the one-loop expression for the bosonic fluctuations re- 
duces aM but is plagued by infrared divergences similar 
to the ones encountered above in the superfluid phase. A 
resolution of this issue is expected from the consideration 
for the momentum dependence of the couplings or from 
the use of functional renormalization group equations. 

As a last application of our formalism we compute the 
speed of sound as extracted from the low energy disper- 
sion relation lu = v m \q\ 18}. In the dimensionless formu- 
lation v s = 2Mv s /kp it takes the form 



v 8 = ^2\ <j) p A <j) . 



(58) 



For the broad resonance limit the result is plotted at 
T = throughout the crossover in fig. |H1 



VIII. LIMITS OF ENHANCED UNIVERSALITY 

In this section we discuss in more detail the status of 
the different limits of enhanced universality presented in 
sect. HV| 
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FIG. 8: Dimensionless speed of sound v 3 for the molecules 
at T = 0. In the extreme BEC regime it approaches v s — > 
2y/c/ (37r). This limit (plotted dashed) is obtained by neglect- 
ing the contribution from the boson density in the equation 
of state. 



A. Exact Narrow Resonance Limit 

A nontrivial exact limit exists for which h<f, — ► and 
\ji —* while c and T are kept fixed. It applies to the 
symmetric phase including the location of the critical 
line. This exact limit remains valid for arbitrary con- 
centration c, even if the scattering length a is arbitrarily 
large. As we have discussed before, the molecules and 
fermionic atoms decouple in the limit — > such that 
the Gaussian (one loop) approximation becomes exact. 
Nevertheless, our limit can describe the full BCS-BEC 
crossover as visible from fig. ^ I n practice, the appli- 
cability of this limit corresponds to a narrow resonance, 
much narrower than the ones currently investigated for 
lithium or potassium. 

We learn from eq. I|17l) that our limit corresponds to 
v — 2a cx hi — ► 0. Due to the simultaneous vanishing 
of v — 2a and h<f, our limit is not the trivial limit of a 
mixture of a noninteracting atom gas plus a noninteract- 
ing molecule gas. Nevertheless, we will see that in the 
narrow resonance limit an appropriate mean field theory 
becomes exact and can describe the high temperature 
phase including the approach to the critical temperature 
of the phase transition. (This mean field theory differs 
from standard BCS mean field theory results, cf. fig. Q]) 
The existence of this limit guarantees that mean field 
theory remains valid as long as remains small, say 
hcf, < 1. This is confirmed in fig. [21 

In order to establish the exact results for this limit 
we first perform in eq. the functional integral for 
the fermions ip. For A^, = this can be done exactly. 
The Gaussian integral yields an intermediate action S[4>] 
depending only on 0, with 

Z= ( V4>e~ m - (59) 



Let us discuss the properties of S[<j>] in the limit h,/, — ► 0. 
The effective potential u in S corresponds to the mean 
field potential and we note that u = — f / (%ttc{o)) + 

uf {f, a) depends on the gap parameter f — h 2 ^(j)*(j) but 

not explicitly on (For the explicit mean field formulae 
cf. sect. [V] eq. (|77|) . and app. [5]) The loop correc- 
tion to the inverse boson propagator P<f,{Q) vanishes for 
h$ — > due to the overall factor hi. In consequence, 
the inverse bosonic propagator for <f> is precisely given 
by the classical part V^ l) = 2mnT + q 2 /AM + v - 2a. 
This is the inverse propagator for free bosons with mass 
2M. All terms in S[</>] with more than two powers of 
4> also involve powers of and therefore vanish in our 
limit. In consequence, the functional integral for <f> be- 
comes Gaussian in the limit Ka = and we can solve the 
functional integral l|59|l exactly. In our limit it should 
be evaluated for <j = D/2. In this limit the number den- 
sity of microscopic molecules Tim therefore corresponds 
to the density of free nonrelativistic bosons of mass 2M 
with vanishing chemical potential l|48[) . as given by the 
canonical partition function. 

In the narrow resonance limit the mean field approx- 
imation for the effective potential becomes exact. This 
yields exact estimates for the critical temperature and 
all quantities of the symmetric phase for T > T c . In 
contrast, the low temperature phase has only a partially 
meaningful limit for — > 0. For nonzero f the bare con- 
densate fraction Qc — 67r 2 f//i^ (cf. sect. IVI|I diverges. 
The constraint tic < 1 then implies that f has to vanish 
cx h/^. In summary, the limit — > corresponds to a 
universal phase diagram in the narrow resonance limit, 
as given by the lowest curve in fig. ^For small h<f, this 
universal curve is smoothly approached (fig. |2J). All cor- 
rections to the "narrow resonance universality" are pro- 
portional to the dimensionless quantity h 2 ^. Obviously, 

for h,f, —>■ all exact results depend only on the param- 
eters c and T, whereby every point in the (c, T)-plane 
may be realized by different combinations of a (or v), n 
(or a) and T. 

At this point a comment on the interpretation of the 
narrow resonance limit is in order. The choice of the clas- 
sical values for the coefficients in the boson propagator on 
the microscopic level corresponds to an intrinsically non- 
local situation. Their actual values are found by simple 
physical considerations, i.e. — 1/4M as the gradient 
coefficient for bosons of mass 2M and chemical potential 
—2a for particle number two. However, for realistic nar- 
row Feshbach resonances additional physical features, as 
a nonvanishing A^,, may become important. 

B. Broad Resonance limit 

The broad resonance limit obtains for h<f, — > oo while 
keeping c fixed. We will see that it corresponds to a 
model for fermionic atoms with local interaction and 
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without explicit molecule degrees of freedom. For — > 
oo all quantities depend only on c and T . The broad res- 
onance limit therefore shows a particularly high degree 
of universality. For fixed c and finite a the limit —> oo 
is accompanied by v — > oo according to eq. (|17|) 



This relation is independent of a. For finite a the in- 
medium scattering length a coincides with the reduced 
scattering length cir (I18fl . Going back to the original 
definition of the action for the functional integral © the 
broad resonance limit corresponds to the double limit 
\v\ — ► oo, — ► oo. In this limit, the "classical" terms 
V0*V</> and <j)*d T (j) become subdominant and can be ne- 
glected. In the fcrmionic language (0) our model there- 
fore reduces to a purely local four-fermion interaction. 
Our partially bosonized description has to match the 
purely fermionic description with a pointlike interaction 
term. In other words, all physical results can only depend 
on the effective coupling (cf. 

- - hi 

\ = \ 4 ,-JL. (61) 

Any individual dependence on h^, D\ or (at fixed A) 
can only be an artefact of an insufficient approximation 
and can therefore be used to test the validity of various 
approximation or truncation schemes. 

As the scale set by h<p drops out in this limit, we can 
express the crossover in terms of a single parameter, the 
inverse dimensionless scattering length c _1 . In the su- 
perfluid phase we can use for the order parameter the 
squared dimensionless fermionic mass gap f = h^ip*^. 
These are precisely the generic parameters used in a 
purely fermionic description. They characterize the sys- 
tem uniquely in the strict limit — > oo, but also remain 
very efficient for hj, 3> 1. Corrections will be 0(h^ 2 ) or 
less. 

The systems currently investigated experimentally are 
6 Li and 40 K. Both range in the broad resonance regime, 
as can be seen from fig. and discussed in sect. [X] 
This clearly motivates the theoretical investigation of this 
regime. 

In contrast to the narrow resonance limit (h^ — > 0, 
A,/, — > 0) the broad resonance limit cannot be solved 
exactly. It still corresponds to an interacting fermion 
mo del, as discussed in detail by Strinati et al. 0, H3, 
Ell 122. l23l| . Nevertheless, it allows us to make a de- 
tailed matching with a purely fermionic description and 
to compare all physical results directly by computing the 
concentration c in both approaches. Furthermore, the 
broad resonance limit can be used as a starting point for 
a systematic investigation of corrections beyond a local 
four-fermion interaction. 

Our present calculations in the broad resonance limit 
still involve quantitative uncertainties related to our ap- 



proximation scheme. This differs from the exactly solv- 
able narrow resonance limit. On the other hand, due to 
the "loss of memory" concerning the details of the micro- 
scopic interaction vertex J3J), this limit has a higher de- 
gree of robustness as compared to the narrow resonance 
limit, where the details of the microscopic interactions 
must be known. 

Finally, we would like to point out that a further 
"crossover problem" , i.e. the crossover from small to 
large h$, emerges from the above discussion. We leave 
this exciting field, physically describing the crossover 
from nonlocal to pointlikc interactions, for future work. 
Our interpretation of the narrow and broad resonance 
limits is confirmed in the numerical study [37| compar- 
ing different types of interaction potentials. 



C. BCS and BEC regimes 

The BCS and BEC regimes correspond to the limits 
c — > 0_ and c — ► 0+. In the limit of small |c|, the 
marginal parameter hj, turns irrelevant. Formally, this 
can be seen form the microscopic four-fermion vertex J3J): 
In these limits, the absolute value of \v — 2a\ becomes 
large compared to the momentum dependent terms. In- 
deed, the temperature and momentum dependent terms 
in the denominator of eq. are typically of the order 
e F = k 2 F /{2M) whereas \V - 2a\ ~ e F h\l\c\. For |c| -> 
the interaction becomes therefore effectively pointlike for 
any given nonzero value of hi. In consequence, the re- 
sults can only depend on A or c, but not on h$ separately. 

Though an intuitively expected result, this kind of uni- 
versality is interesting from a conceptual point of view. 
It is associated to a loss of memory concerning field de- 
grees of freedom: On the BCS side, the molecules could 
have been omitted, and the BCS picture of a weakly in- 
teracting Fermi gas exhibiting the formation of Cooper 
pairs becomes valid. On the BEC side, we end up with 
an effective theory for the "dressed" molecules. They 
behave exactly like weakly interacting pointlike funda- 
mental bosons (more details will be presented in the next 
section), whereas atom degrees of freedom are completely 
negligible for the macroscopic features. 

The reason for this behavior can be seen directly from 
the functional integral J2J) and the microscopic action : 
Introducing the variables c _1 and ip — h^tfi ((<p*)(0) = 
r), which is an exhaustive set to describe these limits, 
we observe the following: For c _1 — > —00 (BCS limit), 
the classical contribution to du/df (|38f) acts as a large 
positive mass term oc — c _1 for the bosons, and the cor- 
relation functions associated to (p are heavily suppressed. 
On the other hand, for cr 1 — > +00 (BEC limit), the so- 
lution of the crossover problem implies for the effective 
chemical potential a — ► — 00, and — cr constitutes a mass 
term for the fermions, similarly suppressing their prop- 
agation. Though the interaction between fermions and 
bosons remains strong through (p^eip^ , the overall role 
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of one of the two different degrees of freedom becomes 
unimportant in the two respective limits. 

It is instructive to discuss the BEC limit from a some- 
what different formal perspective. Let us consider the 
model defined by Sb (|14f> and take the limit a — > — oo. 
The negative effective chemical potential acts as an in- 



frared cutoff 



-a = k , both for the fcrmion and boson 



propagator. We may now integrate out the fermions. 
In the limit k — > oo the contribution from the fermion 
fluctuations to the effective action §[</>] (similar to eq. 
(|59|l discussed for the narrow resonance limit) vanishes. 
Thus S reduces to the action for free bosonic molecules, 

S = /. 4>*(d T - A/2 + v a - 2a)<}). This free theory is 
easily solved and exhibits the standard BEC for T < T c . 

More in detail, we observe for (|27[l that both 70 
and 7 diverge for k — > 00, resulting in 



rf 3 



q 



(2tt 



is 17* " 7 



ATq 2 



(62) 



(2tt) 3 l 2 v g 2 + fc 2 q 2 
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8tt' 



For large k only the remnant of the additive renormaliza- 
tion of v remains 8 . Similarly, for Z§ and A§ one obtains 
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(63) 



Both and A$ take their classical values for k — > 00. 

In consequence, the limit — a — ► 00 is exactly solved. 
Physically, this describes, however, only the very extreme 
BEC limit. Still, for large but finite — a we recover an 
effective bosonic theory that can be solved by by pertur- 
bative methods. 



D. Scaling limit 

A further interesting form of universality was pointed 
out by Ho in the limit c" 1 — > 9 . It is argued that 
if the scattering length drops out, then the density (or 
kp) remains the only relevant scale (besides T), and the 
thermodynamics of the system should be governed by 
simple scaling laws. 

In this form the argument is not complete since it as- 
sumes in addition a strictly pointlike interaction. De- 
viations from the pointlike structure, as reflected in the 



7 The infrared cutoff k should not be confused with the Fermi 
momentum kp. 

8 In the formal limit k — » 00 at fixed momentum cutoff A no 
renormalization of v occurs. 

9 A similar argument was proposed even earlier in the context of 
bosonic systems in l56l . 



Feshbach coupling h<f,, introduce new scales and there- 
fore new dimensionless parameters as h^. This is clearly 
seen by the dependence of the quantities T c and A^ 
(effective gradient coefficient for the dressed bosons) in 
fig. |2 which indeed refers to the resonance c _1 = 0. 
For example, a given corresponds to a particular 

density kp T ^ for the crossover from a broad to a nar- 
row resonance. (We may define kp T ^ by the condition 

h^(h^,,kp r ^) = 10.) In this sense the scaling limit does 
not exhibit complete universality since it depends on an 
additional parameter. 

Nevertheless, for a broad resonance the actual value 
of h<f, is irrelevant for the dressed quantities (as long as 

we are in a range of kp sufficiently away from kp r ^). 
Thus Ho's argument becomes valid in the double limit 
10 c" 1 — > 0, hef, — > 00: For dressed quantities all dimen- 
sionless numbers and ratios can be predicted completely 
independently of the microphysical details of the system! 

Our results for a broad Feshbach resonance extend this 
argument also away from the location of the resonance. 
Universality holds not only for one particular value of B 
at the Feshbach resonance. For arbitrary B in the whole 
crossover region a single parameter c _1 describes all rel- 
evant macroscopic properties of the dressed quantities. 
The concrete microphysics of a system is only needed to 
relate c to B, i.e. it is only reflected in one function c(B) 
which varies from one system to another. 



IX. LOW DENSITY LIMITS 

So far we have developed a rather complete picture 
of the phase diagram for ultracold fermionic atoms in 
terms of the parameters c and h^. In order to make 
contact with experiment we should relate these param- 
eters to observable quantities. In particular, we will see 
how the concentration c is related to the magnetic field. 
One of the great goals for ultracold atom gases is the 
realization of systems with well controlled microscopic 
parameters. This is a central advantage as compared to 
solids or liquids for which usually the precise microscopic 
physics is only poorly known. In order to achieve this 
goal the "microphysical parameters" should be fixed by 
the properties of the individual atoms and molecules, e.g. 
by scattering in the vacuum or by the determination of 
binding energies. The microscopic parameters are then a 
matter of atomic or molecular physics and do not involve 
collective effects of many atoms or molecules. 

Our functional integral approach relates the micro- 
scopic parameters to "macroscopic observables" . It can 
be used for arbitrary values of the density and temper- 



10 In our approach this also holds for the narrow resonance limit 
c _1 — ► O,/i0 — > 0, Aj, — > 0. However, in practice or similar 
terms will not vanish and induce scaling violations. 
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ature. In particular, it can be employed for the compu- 
tation of properties of excitations in the vacuum at zero 
density and temperature. Within one formalism we can 
therefore not only compute the properties of many-body 
systems, but also the scattering behavior of individual 
atoms or molecules in the vacuum. This enables us to re- 
late our microscopic parameters (appearing in the action) 
to observables in atomic or molecular physics. What is 
needed for this purpose is a computation of binding en- 
ergies and scattering cross sections in the limit T — > 0, 
ri — 0. This will be done in the present and next section. 

We are interested in the limit where the density goes 
to zero while T = T/ep is kept fixed at some value 
T > T c . Of course, then also T goes to zero, but we 
do not have to bother with condensation phenomena or 
critical phenomena (for T sufficiently larger than T c ). In 
the limit n — > 0, kp — > also the concentration c van- 
ishes. Actually, we have to consider two limits \c\ — > 
separately, one for positive and the other for negative 
c. For positive c the low density limit corresponds to 
a gas of molecules, while for negative c one obtains a 
gas of non-interacting fermionic atoms. If the scattering 
length is much smaller than the average distance between 
two atoms or molecules, the properties of the dilute gas 
are directly related to the physics of individual atoms or 
molecules. 



A. Atom and molecule phase 

The "fermion gas limit" or "atom phase" (negative c) 
is realized for positive v or negative scattering length 
clr. In this case the molecular binding energy is pos- 
itive such that the excitation of "molecule resonances" 
requires energy. Since for fixed T and kp — > the tem- 
perature becomes very low, T = Tep, the excited states 
become strongly suppressed and can be neglected. In 
this limit one finds a — ► 1 and therefore a — ► ep — ► 
and a — > an. At fixed v and the concentration van- 
ishes c = apikp — > 0. Therefore the interaction effects 
are small. In particular, T c vanishes in the limit c — > 0, 
cf. fig. ^ Finally, the wave function renormalization Z$ 
diverges oc h\J y/jp for a vanishing momentum q 2 — > 
in the propagator. 

For a negative molecular binding energy or v < 
the low density and low temperature state becomes a 
molecule gas. Now the number density of the fermionic 
atoms is suppressed since for zero momentum their en- 
ergy is higher as compared to the molecules. We will 
show that this "molecule phase" is characterized by a 
nonzero negative value of a. 

Let us consider the quadratic term for <j) in the effective 
potential or the term linear in p = <p*<j), 
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According to eq. l|B4ll or 
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the "mass term" 



unless it vanishes. Our limit corresponds to Um cx kp — > 
0, T cx Tkp — > or i%m cx T 3 / 2 -► 0. This requires that 

In consequence, the chemical 



must vanish oc 



kp ■ 



potential a reaches a nonzero negative value 
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dU] 



-(p = 0) and 



- (p = 0) have to be evaluated 



Here — ; — — ; — 

Op ^ ' op 

for a = a a- We arrive at the important conclusion that 
for v < the vacuum or zero density limit does not 
correspond to a vanishing effective chemical potential but 
rather to negative a a < 0. 

The effective scattering length (|17f) becomes in the low 
density limit 
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(67) 



We will see below that (dU\/dp)(p — 0) does not vanish 
for a a < 0, n — > 0, T — > 0. Therefore away from the 
resonance o,{<ja) remains finite for kp — ► and the con- 
centration parameter c vanishes in the low density limit. 
This is consistent with a noninteracting gas of molecules. 
It is not surprising that one encounters the Bose-Einstein 
condensation of free bosons as the temperature is lowered 
below T c . 

In the molecule phase the nonzero a a < adapts the 
additive constant in the energy such that the molecular 
energy level is a zero. In consequence, — a a appears as a 
positive energy (gap or mass term) in the propagator of 
the fermionic atoms. As it should be this suppresses the 
relative number of fermionic atoms in the limit T — > 
and we end with a dilute gas of molecules. The two low 
density limits (with v ^ 0) therefore both correspond to 
the limit of vanishing concentration |c| — > 0. Positive 
and negative c correspond to pure gases of molecules or 
fermionic atoms, respectively. 



B. Binding energy 

In the molecule phase a a has a simple interpretation 
in terms of the molecular binding energy in vacuum, cm, 



11 Eq. 1651 becomes exact if the frequency dependence o f th e cor- 
rection to can be neglected (cf. the appendix in llSl ) and 
+ A t f > (q 2 )q 2 is used to parameterize the exact inverse propa- 
gator at zero frequency. 
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namely 



and therefore 12 for T — > 0, a < 



£m = 2(7,4 • 



(68) 



This can be seen by the computation of the scatter- 
ing amplitudes in app. [D] Alternatively, we may find 
this relation by comparing the densities of molecules 
and fermionic atoms in the low density limit. For small 
enough density and temperature the interactions become 
unimportant and the ratio rip /um should only depend on 
the binding energy in vacuum. Indeed, in our limit a is 
negative and diverges 



lim a 

fcf—>0 



2M<j a 
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(69) 



This simplifies the momentum integrals in the loops 
considerably. The wave function renormalization in the 
"molecule phase" diverges for cm — > 



Z A = 1 
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Furthermore, for large negative a one finds — 
A^/Ztp = 1/2 18]. One therefore has A^/Z^ = 1/4M 
and 7%m corresponds to a gas of dressed molecules with 
mass 2M. On the other hand, the small fraction of 
fermionic atoms obeys eq. (IB4(l with <j = a a resulting in 



^» = ^[exp{I(|l-^} + l] 



(71) 



In the low density limit the interactions become negli- 
gible and np,o should reduce to the density of the free 
gas of atoms with mass M. However, even for zero mo- 
mentum the energy of a single atom does not vanish - it 
is rather given by half the binding energy that is neces- 
sary for the dissociation of a molecule. This corresponds 
precisely to eq. (|71f) . provided we identify ejf = 2a a- 
We conclude that the binding energy is not simply given 
by v (which is defined for a = and in the absence of 
molecule fluctuations) but rather obeys 



£M = ^+ ^p(p = °) G = a A,T = 0). 



(72) 



A measurement of the binding energy of molecules in 
vacuum, ejvf, can relate v to observation. 

For this purpose we investigate the behavior of rfA — 
dU/dp(p — 0) as T — ► 0. First we note that the con- 
tribution from the molecule fluctuations vanishes, since 
for m 2 > one finds (du[ B) / dp){p = 0) oc n M cx T 3 ' 2 . 
For the contribution from the fermionic fluctuations one 
infers from eq. I|38f) (now in dimensionful units) 
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For a < the non-analytic behavior cx y ' —a has an 
interesting consequence for the threshold behavior near 
ejv/ = 0. The equation determining a a in the molecule 
phase (i.e. m|((T^) = 0) reads 

h 2 

v + -l(2Mf/ 2 (-a A ) 1/2 - 2a A - 0. (75) 

Independently of the value of a a and eq. (|75|l can be 
transformed into 



e M — 2a a = 
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We recover the well known universal result for the scat- 
tering length a(aA) [33l l3^ | provided a can be identified 
with the physical scattering length a. This result is in- 
dependent of the value of the Feshbach coupling. One 
concludes tuiy —> 0) — > 0. Close to threshold (a A — > 0), 
the last term in eq. I|75|l becomes subdominant, imply- 
ing a quadratic behavior ey[ oc v 2 . In contrast, far away 
from threshold the binding energy equals the mean field 
binding energy, 6m = v. For broad resonances, the lin- 
ear regime will only be approached far away from the 
resonance when v^A> (2M) 3 / 2 /i 2 /(167r). 



In the "atom phase" we define (cf. app. Q 



em = = )/ Z <t> = v l z <l>- 



(77) 



In our approximation one finds in the atom phase = 1. 
Since cm vanishes for v = we find for the transition 
from the molecule to the atom phase that cm is a con- 
tinuous function of v and therefore of B. The derivative 
dcM/dv is discontinuous at the transition point while 
dajdv remains continuous. 
For both phases, we may write 



2a. 



(78) 



For negative cm the ground state value a — aA is de- 
termined by mMa a) = 0. In contrast, the other limit 
of a gas of fermionic atoms corresponds to a = 0, 
e M = m 2 {a = 0)/Z4, > 0. 

The value of £m could, in principle, differ between the 
two limits n — > with negative or positive cm, due to 
the different ground state values of a. It is interesting 



12 For a > the pole on the real positive axis leads to an additional 
factor i in front of the r.h.s. of eq. 1741 . This has, however, no 
physical impact since a = for the physical solution on the high 
field side of the resonance B > Bg. 
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to address in more detail the question of continuity as 
the magnetic field switches between the two situations. 
In fact, we may view this qualitative change as a "phase 
transition" in vacuum as a function of B. The corre- 
sponding "order parameter" is the value of a at T = 0. 
For B > Bq the "atom phase" is characterized by a = 
and the mass term m| = dU/dp(p = 0, a — 0) is posi- 
tive, corresponding to cm > 0. The "molecule phase" for 
B < Bo shows a nonvanishing order parameter a = a a- 
Now the mass term vanishes, = 0, whereas the atoms 
experience a type of gap \<ta\- The binding energy is neg- 
ative, ejvf < 0. We find a continuous ("second order") 
transition where —a a approaches zero as B approaches 
Bq. In this case one has em (Bo) = and m 2 (£> ) = 
for B approaching Bq either from above or below. (As a 
logical alternative, a discontinuous first order transition 
would correspond to a discontinuous jump of a at Bo- 
Then also £m may jump and not be equal to zero at -Bo.) 
As common for second order phase transitions we find a 
non-analytic behavior at the transition point B = Bo- 



X. PARAMETERS AS A FUNCTION OF 
MAGNETIC FIELD 



So far we have described the universal crossover physics 
in terms of two parameters c and h$. For comparison 
with experiment these parameters have to be related to 
to the microphysical properties of a given atomic system 
and to the external magnetic field B. More precisely, the 
three parameters in the microscopic action QJ, z>a, 
and A^a have to be known as functions of B and the 
effective UV cutoff A. The dependence on A can be elim- 
inated by trading i>\ and h^ t \ for appropriately renor- 
malized quantities v(B) and h^p connected to the atomic 
physics of individual atoms, i.e. evaluated for vanishing 
density and temperature. Furthermore, v can be replaced 
by the reduced scattering length 



a R (B) = - 



4ttD(B) 



(79) 



The residual microscopic pointlike fermion interaction 
\$,A may be related to the background scattering length 

dbg 



MX 



dbg 



4-7T 



(80) 



We will neglect here a possible (weak) dependence of ab g 
on the magnetic field. Then the magnetic field enters 
through olr(B). The microscopic information is now en- 
coded in aji(B), ab g and h^p. We recall that in the broad 
resonance limit the background scattering length db g can 
also be incorporated into ft. 2 (B) . By a suitable Fierz 
transformation we may equivalently use a formulation 
with A^a = 0. 



We need to relate these parameters to our universal 
dimensionless parameters 

c(P,h^,,k F ,a) — * c(B,k F ,cr), 
h<j,,o(i>,h<i„kF,(T) — > h^o(B,k F ,a). (81) 

Here c and h^fi are evaluated in the limit n — > 0, T — > 0. 
Furthermore, the effective Yukawa coupling appearing in 
the relation between the gap f and the squared field p 
includes renormalization effects, /u = h^(h^,fi,a,T,c). 
This is discussed in the appendix 1X1 



A. Concentration 

We may use the measurements of cm(B) in order to 
gain information about our parameters. Consider 6 Li in a 
setting where the "open channel" consists of two atoms 
in the lowest energy states with nuclear spin mj = 1 
and mi = 0, respectively. (These two lowest hyperfine 
states correspond to our two component fermion The 
binding energy near threshold has been measured |57| for 
four different values of B and fits well in this range with 



e M = -P(Bo^B) 2 , 



(82) 



where, for Li, 



/3 (Li) = (7.22 • 10 13 G 2 /cV) _1 = (27.6keV 3 )~\ 



B, 



(Li) _ 



834.1G = 16.29 eV 2 . 



(83) 



The corresponding values for 40 K are 

/3( K ) = (4.97 • 10 9 G 2 /eV) _1 = (1.89 • lO^keV 3 )" 1 , 

(84) 



B { K) = 202. 10G = 3.95 eV 2 . 



From eq. 176|) we infer o,(<ja) = (—^mM)' 1 ^ 2 and we can 
also extract (3 from scattering experiments as described 
in app. El Close to the resonance we can therefore re- 
late the concentration parameter c = a(<7A)kF with the 
magnetic field 13 by 



c-\B) 



fQM\ V 2 

■(^5-) (B-B )=-tb(B-B ). (85) 



The relation (|85|) remains actually also valid for the 
"atom gas", B > B , £m > 0. In this case one has a = 
and 



implies 



D = p(B-Bo) = -^- 



k F c-\B) = Jp^-(B -B). 



(86) 



(87) 



More precisely, c(B) stands for c(B, hp, a a)- 
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It is experimentally verified that the relation a -1 = 
yffJM (B - B) holds with the same value of (3 (jH2l as 
extracted from em and we will confirm this by a "com- 
putation from microphysics" in app. [D] This establishes 
the relation 



16tt 2 /2 2 



such that eqs. lf%5|l and coincide. 

The relation (|85|l permits direct experimental control 
of our concentration parameter for T = and n = 0. For 
6 Li and 40 K one finds in the (arbitrary) unit 14 k F = leV 



= 0.0088G" 1 
t£ = 2.745G" 1 = 



= (113.1G)" 1 
(0.36G)" 1 = 



= 0.45eV" 
140.55eV" 



(89) 



This will be easily extended to T ^ and arbitrary n 
below. The microscopic relation that is independent of 
T and n involves (3 according to 



4:TTV 



((3M) 1 / 2 (B-B ) = T B k F (B-B ). (90) 



We emphasize that the determination of t b in eq. I|85|) 
does not require knowledge of since we can use directly 
the experimental determination of j3 from eq. 1)82(1 . 

Away from the resonance the issue of the correct choice 
of c is more subtle. The observed scattering length a is 
well approximated by a resonance term 



1 



1 



B 



(91) 



and a constant (or weakly B 
scattering db g , 



dependent) background 



a{B) = a res (B) + a bg . 



(92) 



The question arises how the background scattering is 
properly incorporated. In principle, not only the ex- 
change of the dressed molecules contributes to the scat- 
tering, but also a residual atom interaction in vacuum 
which is parameterized by the effective coupling A^o- I n 
the broad resonance limit our results are consistent with 
a minor role of A^o in the molecule phase and we will 
neglect it. 

Indeed, we will see below how the background scat- 
tering is incorporated into cl{cta) for the molecule phase, 
such that we can associate cl(<ja) with the observed scat- 
tering length, including the background term 



a(a A ) = a(B). 



(93) 



For the concentration this implies a simple relation to 
a and therefore a somewhat more involved dependence 



on B. Our central relation between the concentration c 
and the magnetic field is simply expressed in terms of the 
measured scattering length a(B), 



a(B)ki 



(94) 



This relation is quite satisfactory. Indeed, in the broad 
resonance limit the interaction is essentially pointlikc, 
both for abg and a res . Physical results can therefore only 
depend on a = a res +ab g - This holds independently of the 
precise value of A^a- Indeed, by a Fierz transformation 
we can shift the interaction from A,/,, a to the contribu- 
tion from molecule exchange ~ h^/v and vice versa [Isj . 
Since for pointlike interactions the Fierz transformations 
are exact identities the physical results should only de- 
pend on the sum A / — 2a a) — A^a and not on the 
value of an individual piece like A^a- In terms of the 
renormalized quantities the results can then only depend 
on a and not on a res and at g separately. This feature 
is precisely realized by the relation 194(1 : all results de- 
pend only on the concentration c, and this is in turn is 
given by the total scattering length a, independent of a 
decomposition into a res and a^g- 

Based in the above argument we will use the simple re- 
lation 1(94(1 also for B > Bq- Within our approximations 
this implies that we associate an = a res + ab g for B > Bo 
which amounts to B - dependent h^ (B). In fact, our 
approximation does not resolve the Fierz ambiguity for 
B > Bo- In app. [D]wc discuss in detail the scattering of 
atoms in vacuum. There we comment on the issue of the 
Fierz ambiguity more extensively. 

In summary, the dependence of c _1 on B obeys the 
nonlinear relation 



c 



1 

d res 



1 



1 -\- CLfjg I &res 

t b (B-B ) 



(95) 



1 - a bg k F T B (B - Bo) 



Of course, kp is not related to n in this case. 



where f B = k F r B = y/JJM = 0.45(140.55)eV for 6 Li 
( 40 K) is independent of k F - This nonlinearity in the re- 
lation between the magnetic field B — Bo and the cou- 
pling strength c~ 4 is an important effect. For 6 Li this 
produces the strong increase of the bare molecule frac- 
tion on the low field side of the resonance as seen in fig. 
13 We note that for negative ab g (as for 6 Li) the scatter- 
ing length reaches zero at some value of B in the BEC 
regime, B - B = -A = aj^jSM) -1 / 2 . At this point c 
vanishes. For B below Bo — A the system belongs to the 
weakly coupled BCS regime. The transition between the 
BCS and BEC regimes in the weak coupling limit can 
therefore be realized by experiment. For 6 Li one finds 
A = 300G, for 40 K one has A = -7.8G. 



B. Yukawa coupling 

The discussion of the relevant range for our second pa- 
rameter, the dimensionless Yukawa coupling h$, is more 
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involved. Let us denote the effective Yukawa coupling 
in vacuum in the atom phase by Ln or h^y for the di- 
mensionless counterpart. In app. 0we have computed 
its relation to the microscopic coupling /i^a appearing 
in the action QJ. From eq. (|88|1 and the experimental 
determination of [3 in the atom phase we can extract the 
ratio 



4tt 



/x (M 3 /3)V2 

= 2.516 •l(T 7 eV~ 2 G = 4.915 • 1CT 9 . 



(96) 



(Here the numerical value is given for 6 Li.) For the broad 
Feshbach resonance in the 6 Li system the molecule state 
belongs to a singlet of the electron spin, resulting in [Im ss 
0. The microscopic value of p, = di>\/dB obtains its 
essential contribution from the magnetic moment of the 
atoms in the open channel, which is well approximated 
[58| by the Bohr magneton /i g = 5.788 • 10- n MeV/T = 
0.2963MeV _1 . For 40 K we take ft = 1.57 fx B HJ. With 

/i (Li) = 2/x B , fiW = 1.57/i B , (97) 

and kp — leV this yields high values 

610, (h% ) 2 = 3.72 ■ 10 5 , 

(h*) 2 = 6.1 • 10 3 . (98) 



h u 

hlo = 79, 



Using the broad resonance criterion h 2 ^ > 100, we see 
that indeed 6 Li and 40 K belong to this class, cf. fig. 

In the molecule phase the vacuum value of differs 
from h^,fi, as we have shown in app. ^ If we incorporate 
dbg in A^, ; a such that /i^a is independent of B one obtains 



0,0 ~~ 



10,0 



1 - a bg /a((JA) 



(99) 



It is interesting to see how the renormalization effects 
turn <i((ja) into the full scattering length including the 
background scattering. Indeed, assuming that h^fl is in- 
dependent of B and v oc B — Bo we may identify the 
resonant contribution to the scattering length as 



4ttis 

For a(a A ) this yields 



B-Bn 



(100) 



a(a A ) 



m to; 

47r v — 2a a 



ares(^) .(101) 



Precisely at the resonance the renormalization effects 
vanish. 

Away from the resonance, however, we have to take 
into account two renormalization effects which appear in 
the order ab g /a res . First, the ratio (h ™ /h^.o) 2 yields a 
B dependent factor (1 — a^g/a) 2 (with a = a{aA))- Sec- 
ond, the "renormalization of the detuning" by the shift 



v — > v — 2a a induces an additional B - dependence by a 
factor v j[p — 2a a) = a res /a. Eq. I|101f> yields 



a(l - a bg /a,y 



/ d res \ 

H- 

\Cl- dbgS 



(102) 



and we conclude that a is indeed given by a res {B) + a^g. 
We observe that a bg is related here to the microscopic 
"background atom interaction" A^, via the renormaliza- 
tion of h^fl discussed in the appendix lAl Thus eq. I|102|) 
gives a microphysical explanation for the nonlinear B - 
dependence of a. 

This issue is different for the atom phase where the 
B - dependent renormalization effects for h^fl and v are 
absent. In this case a = an(B) only reflects the resonant 
contribution to the scattering if we keep indepen- 
dent of B. The background contribution has then to be 
accounted for by the direct coupling A^.o- Alternatively, 
we may absorb A^a into the molecule exchange contribu- 
tion by a Fierz transformation, including an effective B 
- dependence of h 2 ^ for B > Bq. Then an has to reflect 
the full dependence on B, including the background scat- 
tering. If A^,o can be neglected this requires clr — a(B) 
and leads to the relation (|94() . 

The determination of ejvf = 2aA amounts to the con- 
dition 

™%M v~2a A dUt 

-™- = -n- + - gF (r = 0,^) = 0. (103) 



The bosonic fluctuation contribution vanishes for T — > 
and dU[ F " > /df only depends on a, but not on h$ or 
p,. This is the reason why the behavior for B — > Bo 
can be used to determine c and the ratio l|96|l . but not 
hef, separately. In principle, independent experimental 
information on h<f, t o can be gained from measurements 
of the binding energy outside the resonance. Indeed, for 
(-M < the general relation between cm and v reads 



7 (1 -y/1- 29/ j), 



(104) 



hjpM 3 
32ir 2 ' 



Measurement of cm{B) over a wide enough range allows 
for the extraction of both h 2 ^ and v independently. Far 
away from the resonance one has \"f/v\ <C 1 such that 
(-M = v. From this range we extract the B - dependence 
of v, dtu /9B — dv/dB = fi. On the other hand, close to 
the resonance the opposite limit \~y/v\ 3> 1 applies. Then 
6m = —v 2 l(2"{) involves the concentration c, and, using 



/i(B—Bo), one can extract h^o from 7 



7(2£m). 



Alternatively, one may extract v(B) and h^^(B) from 
a detailed investigation of the atom scattering near the 
resonance, as discussed in app. [D] 
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General relation between c and B 



In subsect. A we have related the mean field scat- 
tering length an in vacuum with the detuning magnetic 
field B — B . This can be used directly for establishing 
the relation between the concentration c and B — Bq for 
nonzero T and n, using the definition (|17f) 



= k 



hiM 



One finds 



1 16na 



Employing the relation 



-t b (B-B ) 



(105) 



(106) 



(107) 



we may also write eq. (110511 in the form 

2Mfl(B - B ) 



v = 2a- 



K 

8ttc 



(108) 



We observe that the dimensionless Yukawa couplin g hj 
appears in eq. (|106|1 . It is related to /i0 jO by eq. (|A6|) 
in app. ^ For practical calculations we use the relation 
(|108jl between v and B — B . 

In the broad resonance limit for large values of /i| of 
the same order of magnitude as the vacuum values l|98|) 
the terms ~ h^ 2 in eq. I|106|l can be neglected if |ct| is 
not too large. This yields the simple expression 

c- 1 = -b. (109) 

We recall, however, that this formula should not be used 
far in the BEC regime where a takes large negative values 
which diverge in the vacuum limit. For moderate or small 
values of hff, the term oc a in eq. (|106fl should always be 
included. 



XI. DRESSED AND BARE MOLECULES: 
COMPARISON TO EXPERIMENTS 

A. Bare molecules 

The fraction of closed channel molecules is given by 
&B = &m + ^c- It has been measured by a laser 
probe that induces a transition to an excited molecu- 
lar level - the quantity Z in [{| equals fls- Indeed, the 
laser probe couples directly to the total number of "bare 
molecules" {4>*<j)) — tib- Importantly, we have seen in 
sect. Ivl and fvTI that this quantity involves the Yukawa 
coupling explicitly, cf. eq. In fig. |3]we com- 

pare our results with the measured value for Li, using 
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FIG. 9: Contributions Oc, &m to the fraction of closed chan- 
nel molecules Qb — Oc + Om- We plot the result for 6 Li 
[h^fi = 610) for T = 0. The condensed part dominates fls 
in the BEC regime, while in the vicinity of the resonance, the 
contributions are comparable. The dashed line for Oc omits 
for c" 1 > the background contribution to a and reflects the 
Fierz ambiguity. 



T = 0, k F = 0.493eV=250nK. The solid line uses 
c = kF(a res (B) + dbg) H93fl for the whole range of B. 
The agreement with the measured value is very convinc- 
ing. In order to give and idea of the Fierz ambiguity for 
B > Bq we also display by the long-dashed line the re- 
sult of the choice c = kpa res (B). For B < Bq the Fierz 
ambiguity is removed by the renormalization effects as 
discussed in sect. IX Bl We emphasize that the inclusion 
of the renormalization 1A12(1 for the Yukawa coupling is 
crucial. Omitting this effect (dashed line) for B < Bq in 
fig. |3] results in a clear discrepancy from the observations 
sufficiently far away from the resonance. 

At this place we may give a few more details of our 
computation of Qb , 



n, 



n 



M 



3tt 2 (- 



2^- 

hl 
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We use the Yukawa coupling as defined in (|A6|) in the 
place where it appears explicitly. All other quantities are 
taken from the solution of the equations determining the 
crossover problem in the superfluid phase, eqs. I|B9IB11|I . 
The formulae for tic, 



6n* 

F 



|0o| 2 =6Tr 2 (l-a bg /a) 2 j 



(in) 



relates the closed channel condensate fraction to the su- 
perfluid order parameter 4>q. In that sense Partrigde et 
al. indirectly measure the superfluid order parameter! 
From VLc we can infer CIb = ^c(l + ^m/^c) by ex- 
tracting the ratio Hm / from fig. In the BEC limit 
Hm/^c becomes negligible, while at resonance (c^ 1 = 0) 
we find fi M /O c = 1-89. 

Let us briefly comment on the density dependence of 
our result. This is particularly simple in our dimension- 
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less formulation. Since fo is taken from the solution of 
the crossover problem it does not depend on kp if c is 
kept fixed. However, in order to find the right dimen- 

~ — 1 /2 

sionless value h$ = 2Mh c j > /kp , kp must be inserted. 
For fixed c this yields tic oc kp. The ratio ti,M/tic de- 
pends also on c. The dependence of tig on c introduces 
an additional density dependence of tig if the scattering 
length a is kept fixed. (Only for fixed akp one always 
has tig oc kp in the broad resonance limit.) 

In the scaling limit cr 1 = the value fo and the ratio 
Qm/&c are independent of kp if we assume the broad 
resonance limit — > oo. In this limit we thus can con- 
firm the simple scaling law at resonance advocated by 
Ho 01 1 oc kp. A similar result has been obtained by 
Levin et al. , though neglecting the contribution ti,M 
which is of the same order of magnitude at the resonance, 
cf. fig. H 

In the BEC limit S1m/^c becomes negligible. Fur- 
thermore, since f oc c _1 oc kp , one finds that tis be- 
comes independent of kp. Note that this result is only 
valid for kp which respect the BEC regime condition, 
cT 1 = (akp)^ 1 > 1 or kp < a -1 . 

In the deep BCS regime c _1 = |a/cf| _1 ^ 1, we also 
find numerically that Hm/Hc <C 1. Using the stan- 
dard BCS result relating the superfluid order parame- 
ter at T = and c , our scaling form yields tig oc 
kp exp(— 7r/ (akp)). Again, in order to stay in the desired 
regime, kp must be restricted to values kp <C |a| _1 . In 
sum, we find the following scaling behaviors with kp 

{const. BEC 
kp Resonance . (112) 

k F e-*/( ak *-) BCS 



B. Dressed Molecules 
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FIG. 10: Condensate fraction Qc and total fraction of dressed 
molecules Q B = SIm + ^ c at T = 0, for 6 Li (a) and 40 K 
(b). Due to the nonlinear relation of c" 1 and B — Bo which 
is particularly strong for 6 Li due to the large background 
scattering value, the curves are twisted compared to the 
f^c" 1 ) plots. We adjust k F to the values found in yj (&F = 
1.50eV=(2500a B )" 1 ) and Q (k F = 1.39eV=(2700aB) _1 ). 
The dashed line reflects the Fierz ambiguity as in figs. 13191 



The condensate fraction measured in 0, qualita- 
tively refers to the condensation of "dressed molecules" 
or di-atom states. At the present stage, however, the 
precise relation between the measured observables and 
the condensate fraction fie has not yet been established 
with sufficient quantitative accuracy. As argued in sect. 
Ivl flc is given by the expectation value of the renormal- 
ized field 

lie — 2 — — = 67T Za,p — 07T p, (113) 

n 

4>R = z x J 2 ^. 

In fig. ^)| we plot our result for the condensate fraction 
of dressed molecules at zero temperature as a function 
of magnetic field, for both the 6 Li (fig. ^| (a)) and the 
40 K system (fig. ^| (b)). Unlike the condensate frac- 
tion of bare molecules, the contribution from the dressed 
condensate is an 0(1) quantity. 

We find qualitative agreement with the observations 
0, Q whereas for a quantitative comparison one would 



need a more accurate relation between the measured ob- 
servables and the condensate fraction flc as defined in 
our setting. First, the scales of the magnetic field B 
for which flc decreases from rather large values to small 
values match the scales found in the experiments. This 
is another confirmation of our universal relation between 
B, the particle density (or kp) and the quantities c" 1 , h^. 
We find that in the BEC limit, the condensate fraction 
approaches 1 very slowly. The observed condensate de- 
pletion is the effect of the interaction between the dressed 
molecules, as expected for a weakly interacting Bogoli- 
ubov gas. 

An interesting quantity is the value of the condensate 
fraction at the location of the resonance, B = Bq. For 
T = we find a universal value f2£, = tic(T — c~ x = 
0) = 0.30, which therefore should apply both 15 for 6 Li 



15 If one would identify the "condensate fraction" quoted in refs. Q, 
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and 40 K. In order to judge the reliability of this universal 
result we may consider the relation 



n c = 



67T 2 f 



(114) 



which involves the renormalized Yukawa coupling hi = 
h^/Ztp. In the broad resonance limit and for T — the 
value fg = fo(T = c _1 = 0) = 0.28 is universal. The size 
of Vl c is therefore determined by h^, Vl c = 16.85/i^ 2 . In 
the same limit we find from eq. I|29|) 



h 6 2 = 7TJ> I d ^ 



q 2 (q 2 - a) 



[(q 2 - a) 2 + f ] 3/2 



(115) 



where a = a(T = c" 1 = 0) = 0.50. This yields hT 2 = 
0.018, n c = 0.30 as seen in fig. \M Away from the 
resonance both fo and a depend on c. In fig. Illl we show 
/i 2 as a function of c _1 for T = 0. 

In the broad resonance limit the value of the renor- 
malized Yukawa coupling only depends on c _1 and T. 
This universal value hi(c ,T) is reminiscent of the ex- 
istence of a partial infrared fixed point as some suitable 
infrared cutoff scale is lowered. The value of h<j, at the 
fixed point typically depends on c _1 and T, but other- 
wise the memory of the initial conditions of the flow (e.g. 
/i0,a) is lost. It is also plausible that also is governed 
by some fixed point. If this picture is true, the precise 
location of the fixed point may change as the approxi- 
mation method is improved, for example by taking the 
molecule fluctuations into account for the computation 
of and the renormalization effects for ft, 2 . Neverthe- 
less, if both the 6 Li and the 40 K system are within the 
range of attraction of the fixed point, the value of Vic 
for T = c _1 = must be the same! Inversely, if the 
observations should establish different values of £l c this 
would imply that either no such fixed point exists or that 
one of the systems is not yet close enough to the fixed 
point. Only in this case the system could keep additional 
memory of the microscopic Yukawa coupling (e.g. h^.o) 
or the background scattering length ai, g such that Q c 
could depend on these parameters. Our findings so far, 
however, strongly suggest the universal behavior predict- 
ing for broad resonances a universal curve Slc^c -1 ) for 
T = 0. 



XII. CONCLUSIONS 

The theoretical description of the crossover problem 
for ultracold fermionic atoms can be split into two steps. 
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FIG. 11: Renormalized Yukawa coupling hi as a function of 
c _1 at zero temperature in the broad resonance limit. 



First, the computation of universal properties in terms of 
the two relevant dimcnsionlcss couplings c and h$. Sec- 
ond, the relation of c and h^ to the particular parameters 
of a given Feshbach resonance as the detuning of the mag- 
netic field B — B . Both steps can be performed within a 
common setting of non-relativistic quantum field theory, 
described here in terms of a functional integral. 

The first step can be done for all densities and temper- 
atures and for arbitrary Feshbach resonances that can be 
approximated by our system of one open and one closed 
channel. It reveals the universal aspects since all observ- 
ables are expressed in terms of only three parameters 
c, h$, T. For broad Feshbach resonances (large h^) the 
universality is further enhanced. The renormalized quan- 
tities relating to the "macroscopic" properties of dressed 
molecules and atoms only involve a renormalized Yukawa 
coupling /i^ that takes a "fixed point" value h<f,(c, T) that 
is computable in terms of c and T. The "dressed quanti- 
ties" therefore depend only on c and T. In consequence, 
for T = or T = T c they turn to functions depending 
only on the concentration c. Furthermore, in the scaling 
limit at resonance (|c| — > oo) the renormalized observ- 
ables depend only on T and take predictable universal 
values for T = or T = T c . 

As a particular aspect of universality we find perfect 
equivalence between a purely fermionic description of 
the crossover problem and a model that uses explicit 
bosonic degrees of freedom corresponding to the micro- 
scopic molecules in the closed channel of the Feshbach 
problem. The microscopic molecules are a convenient 
tool to take into account the non-locality of the effective 
fermion-interaction induced by the molecule exchange. 
This nonlocality is relevant for narrow and intermediate 
Feshbach resonances. For broad Feshbach resonances, 
however, the nonlocality becomes subleading as far as the 
properties of the dressed atoms and molecules are con- 
cerned. The broad resonance limit h$ — > oo precisely cor- 
responds to a pointlike microscopic fermion-interaction. 
As it should be, only one dimensionless parameter (be- 
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sides T) characterizes the system, namely the scattering 
length a appearing in the dimensionless concentration 
c = akp. 

For a test of our universality hypothesis we have 
computed several macroscopic observables as functions 
of c,hcf,,T. This concerns the fraction of microscopic 
molecules in the atom gas Cls, the condensate fraction 
Qc an d the fraction of uncondensed dressed molecules 
We also compute the correlation length for the 
molecules £m as well as their scattering length in medium 
ajv/. The low momentum excitations in the superfiuid 
phase are described by a Bogoliubov theory for dressed 
molecules even if the microscopic molecules play no role. 
We have computed the corresponding sound velocity vm ■ 
All these quantities may be experimentally tested. 

For these tests we need the second step, the relation 
of c, h<j, and T to the physical properties of experiments 
with a given Feshbach resonance at a given density n and 
temperature T. Here T — Tjtp only involves T and n, 
€p = k F /(2M), n = kp/(3TT 2 ). In order to relate c and 
to observables in the two atom system (at zero den- 
sity and temperature) we compute the molecular binding 
energy cm and the cross section for the atom-atom scat- 
tering a(s) as functions of c and h^. (Here s denotes the 
momentum of the scattering atoms.) This is achieved by 
simply taking in our formalism the limit n — > 0,T — > 0. 
Comparing ejw(c, /i^), cr(s, c, /i<^) to the measured values 
as functions of the magnetic held, cm(B), <t(s, B), we can 
extract c(B) and h^B). In particular, for the broad reso- 
nance limit c(B) is simply related to the scattering length 
a(B) describing a pointlike interaction, c{B) — a(B)kp. 
For a broad Feshbach resonance and "dressed quantities" 
the only experimental input needed is therefore a(B). 
Only for "bare quantities" (as the fraction of microscopic 
molecules) one needs, in addition, the determination of 

The broad resonance limit leads to striking predictions 
that challenge experiment. For example, the condensate 
fraction Qc should take the same value at resonance, 
B = Bq, for all broad Feshbach resonances. (Here Bq 
corresponds to a vanishing binding energy in vacuum, 
cm (Bo) = 0). For T -> or T -> T c this should hold for 
arbitrary density, whereas more generally flc depends on 
Tjtp and therefore on n in a universal way. In our ap- 
proximation we hnd Qc — 0.3 for T = but the precise 
value could change somewhat for improved approxima- 
tions. If both the broad resonances for 6 Li and 40 K are 
well described by the broad resonance limit the experi- 
ments should hnd the same fie*! While the present 
experiments seem to indicate a substantially higher con- 
densate fraction for 6 Li as compared to 40 K, it remains 
to be seen if this tendency remains once the detailed re- 
lation between the observables in both experiments and 
fic is established. 

The broad resonance limit poses the challenge of a 
strongly interacting quantum held theory. Several issues 
wait for an improved treatment. The understanding of 
the renormalization flow for and A^, should reveal the 



existence of an infrared stable (partial) fixed point and 
should allow to judge if both 6 Li and 40 K are within 
the range of attraction of this fixed point. Indeed, the 
universal value of Q.cic 1 = T = 0) requires that for 
both systems the value of h 2 ^ is sufficiently close to the 
fixed point value. The functional renormalization group 
[52l l53l l54j seems to be an appropriate method to ad- 
dress these questions. This should also include the effect 
of molecule fluctuations in the renormalization flow of 
/i^, Z$ and A,/, as well as the computation of the missing 
corrections to the wave function renormalization and gra- 
dient term for the fermionic atoms, and A^,. Another 
issue that is not yet settled in a satisfactory way in the 
present treatment concerns the behavior in the superfiuid 
phase for T — > T c and the universal critical behavior in 
the vicinity of T c . These issues are mainly related to the 
running of the molecule interaction A^. 

Despite the present shortcomings we hope that our 
functional integral approach opens the way for a system- 
atic quantitative understanding of the relation between 
measured microphysical parameters and macroscopic 
collective phenomena for ultracold fermionic atoms. 
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APPENDIX A: RENORMALIZATION OF 
FESHBACH COUPLING 

So far we have not been very precise about the mean- 
ing of the Yukawa coupling h&. Just as v it has to be 
renormalized, this time with multiplicative renormaliza- 
tion. In the action appears the microscopic Yukawa 
coupling ft.0.A (and similar for A^,.a). In the vacuum 
(T = 0, n = 0) one can measure the Yukawa coupling 
by scattering experiments of single atoms. This involves 
a UV - renormalized coupling corresponding to the effec- 
tive molecule - atom vertex at low momentum. We will 
use the vacuum value in the atom phase, h<j,fl, in order 
to define a UV - renormalized Yukawa coupling. Finally, 
in an atom gas at arbitrary T, n the relation between 
the gap and </>o involves a renormalized Yukawa coupling 
h$(T, a). Corresponding dimensionless quantities are de- 
noted h^^ih^fiih^iTju). The quantity h<f,(T, a) typi- 
cally appears in the fermion loop through the gap r and 
we use the shorthand h^. Actually, the situation is even 
more complex since a momentum dependent h ( p(q) ap- 
pears in the momentum integrals of the Schwinger-Dyson 
equations. We will neglect this in a simplified treatment. 

We estimate the UV-renormalization effects for the 
Yukawa or Feshbach coupling by the use of an appro- 
priate Schwinger-Dyson equation shown in fig. 1121 This 
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FIG. 12: Graphical representation ol the Schwinger-Dyson 
equation for the Yukawa or Feshbach coupling. Fermions are 
represented by solid, molecules by dashed lines. The vacuum 
value of the Yukawa coupling, h^fi, is signalled by the small 
blob, the full quantity by the fat shaded blob. The vacuum 
value of the background scattering vertex \^ enters as a small 
square. 



yields 



1 f d 3 q tanh70 
4TJ (2tt 3 ) 70 



(Al) 



Jh.A = — 



We first evaluate Jh.A for T — ► 0, a — > and find that 
the integral is proportional to the ultraviolet cutoff A 



J htA (T = 0,a = 0) = 
Thus the relation 16 between h^ t o 
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and /i<^,a involves A 



(A2) 



hc/,,0 = h<f,,A — h^Xip^A 
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Next we express in terms of h^.o instead of /i^a 
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where Jh converges now in the ultraviolet 
1 f d 3 q rtanh 70 ATM 
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All this is completely analogous to the renormalization 
of &A1& discussed in detail in where the same 

momentum integral appears (up to a prefactor). The 
relevant momentum scale for is now low, adapted 
to the momentum range dominating Jh- We employ 
\<p = Ancibg/M and obtain for the dimensionless coupling 
(a bg = -0.38cV for 6 Li) 



h 



0,0 



1 + 8ira bg k F J h 



(A6) 



with 
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For the symmetric phase a first interesting limit con- 
cerns the BEC regime for a — > — oo where 



lim J h = — 

tJ * oo (57T 



such that 



hd, = hd,o 
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A second limit relevant for the BCS regime takes f — > 
for a = 1, 



lim J h (a 



1) oc InT. 



(A10) 



These two different limits have important consequences 
for the value of the renormalized Yukawa coupling in vac- 
uum. For the atom phase (positive binding energy ej\j) 
we conclude that the limit kp —* 0,T — ► discussed in 
sect. IIXI implies a vanishing of the correction ~ kpJh 
in eq. (|A6|) such that the vacuum value of the Yukawa 
coupling is indeed given by h^o. In contrast, for the 
molecule phase (cm < 0) the limit (|A8|) applies. With 
6m = 2cr a = — 1/(Mo, 2 ((ta)) l|76(l one finds an enhance- 
ment of the vacuum Yukawa coupling for a^g > 
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(m) 
0,0 



l 0,O 
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This correction vanishes at resonance. 

We may also evaluate eq. (|A6ll in the superfluid phase. 
Using the standard BCS gap equation f (|B9|l with VLm = 
0), one obtains 



l 0,O 



CLbgkpC 
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This formula should be reliable for low T where indeed 
fiju small. Eq. (|A12|) smoothly connects to the vacuum, 
eq. (lATTJl . Despite several shortcomings in the derivation 
of i|Allfl . mainly due to the neglect of the momentum 
dependence, we believe that the relation (|A11(I should 
hold approximately also for a more reliable treatment as, 
e.g., the functional renormalization group equations. 



APPENDIX B: GAP EQUATIONS 



1. Normal phase 

Since the loop integral is dominated by high momenta q 2 f» A 2 

we use here X^ A (instead of a more accurate description in terms In the normal or symmetric phase the superfluid gap 
of a momentum dependent A^, (ij 2 )). vanishes, po = 0, and therefore the mass matrix 147|) 
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becomes degenerate. The value of is determined by 
a gap equation (TIm depends on ml, but not on A^), 



m 2 _ m (F)2, 1 \ (F) 0» r 



with 



m 



(F)2 



v + Am 



(F)2 



(Bl) 



(B2) 



evaluated at po = 0. Here we have introduced a renor- 
malized detuning v — v /Z^ and work in the broad reso- 
nance limit such that the term —lojZ^ = Oih^ 2 ) in the 
classical contribution to the boson mass can be neglected. 
Explicitly, we use 
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^r7 3 [tanli7 — 7 cosh 2 7] . 
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The effective chemical potential a is determined by 
the density equation Qf,q + = 1 (fif,o = in the 
symmetric phase), with 
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The polylogarithmic function Life(z) obeys 
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For the computation of SIm we finally need which 
reads in the broad resonance limit 



d 3 q 

48T 3 J (2tt) 3 



q 2 j 3 [ tanh 7 — 7 cosh 2 7] . 
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All displayed expressions are effectively - indepen- 
dent in the broad resonance limit since the renormalized 

l/2~ 

Yukawa coupling htj, — Z^ becomes independent of 

h<f, for hcf, — > 00. Using eq. H48|) we can relate v to the 
concentration c 



v v 



1 
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With the condition flp,o + = 1 we can determine cx 
and m^ self-consistently. 

A second Schwinger-Dyson equation is obtained for the 
self-consistent determination of the dimensionless renor- 
malized four-boson coupling A^ |18|. namely 



A«i = 



A 



1-/a' 



(B8) 



/a = ^/ /^ , Q a -1 [Yexp 2a — l) + 2a sinh -2 al , 

2t7(2tt) 3 lv 7 J 

+ 

a = = -. 

2T 

Approaching the critical temperature from above, 
m^(T — > T c ) — > 0, this leads to an infrared divergence 
of the fluctuation integral / in the second line. In turn, 
this will lead to a behavior A^T — > T c ) cx m^ — * 0. This 
"infrared freedom" at the phase transition is a generic 
feature of bosonic systems. 



2. Superfiuid phase 

The superfiuid phase is characterized by a nonvan- 
ishing superfiuid order parameter po or (squared) gap 
f — h^py. The minimum condition yields a gap equation 
in the broad resonance limit 



Amf> 2 + JLA^ M = 



(B9) 



which expresses the vanishing of the full mass matrix 
term m 2 = (Goldstone theorem) . The fluctuation inte- 
grals entering (|B9|I are the same as in (|B3f) , with the only 
difference that they have to be evaluated at po ^ 0. This 
makes these expressions more complicated. For their ex- 
plicit form, we refer to . Dividing eq. I|B9|) by h 2 ^ we 
find 
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8^c 



2 ( Am 



(F)2 



h 2 



(BIO) 
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(F) f d 3 q A^q 2 + A 0( o o /2 



(2tt) 3 ^<f + 2A 0( o o 



x(exp^A 4> q 2 + 2\^p /T -l) . 

Note that without the bosonic fluctuation contribution 
oc fijif, eq. ijBlOfl is just the usual BCS gap equation 
as obtained from a purely fermionic approach, since the 
factor h^ 2 cancels the corresponding factor in (|B4(1 . Due 

to the additional term cx A^ our approximation scheme 
goes beyond the presently available purely fermionic ap- 
proaches even in the broad resonance limit, since it offers 
the possibility to include fluctuations of the dressed bo- 
son field as well! However, the quantitative treatment of 
the bosonic fluctuations still involves substantial uncer- 
tainties, especially for T near T c in the superfiuid phase. 
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In general, the boson fluctuations are expected to af- 
fect our quantitative results mainly close to the critical 
temperature. At very low T, they should be negligible. 
Therefore, our results should be reliable in this regime, 
which is confirmed by the matching of recent QMC re- 
sults |H3| with our approach [Isj . 
The atom densities obey 



The two lines representing the critical temperature in 
fig. n are obtained by either including (short dashed) 
or omitting (solid) the boson fluctuations in eq. USE}. 
Other current approaches, e.g. [22I ItH lis j . always work 
with the standard gap equation omitting the boson fluc- 
tuations. 
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For a broad resonance the condensate fraction £lc = 
nc/n = 67r 2 po = 67r 2 f//i0 becomes independent of the 
the value of hs- For the actual computations we use 
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where flc can be safely neglected for a broad resonance. 
Since J7m depends on A^ in the superfluid phase, we need 
another gap equation for the four-boson coupling A^. We 
find \$ = for the infrared value in the superfluid phase, 
and we refer to JlSj] for a discussion of how we actually 
deal with this coupling. 



3. Phase boundary 

At the phase boundary, the system is governed by par- 
ticularly simple equations since mi = A^ = po = 0. The 
vanishing of m 2 ^ implies the constraint 

-i4K 2 + i^ M )=0 (B13) 



and the density equation can be expressed as 
Qf,o + £Im = 1 ■ 



(B14) 
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The fermionic integrals for Am^ ' , A^ , AA^ are the 



same as in eq. I|B3|) with Qm simplified to 
3r(3/2)C(3/2) / f \3/2 
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At the phase boundary, the solution of the system 
(|B13IB14|I amounts to the self-consistent determination 
of a, T c . The critical line therefore indeed depends only 
on c in the broad resonance limit hs — > 00. 



APPENDIX C: DISPERSION RELATION FOR 
MOLECULES IN MAGNETIC FIELD 

Due to the interactions between unbound atoms and 
molecules the relation between energy and momentum 
(or frequency and wave vector) for single molecules dif- 
fers from the case of free molecules of mass 2M where 
u> = q 2 /AM. We can extract the dispersion relation Ld(q) 
from our approximation to the full molecule propagator 
'P<p{Q) = 7^)(go, <t) in vacuum, where n = T = 0. In the 
absence of spontaneous symmetry breaking we work in 
the (jf , <p basis where V multiplies the unit matrix. The 
dispersion relation uj(q) follows from 



P^kj,q) = 0. 



(CI) 



For the "molecule phase" (cm < 0) one has to use o = a a 
with cm = 2a a such that 7^,(0,0) = 0, whereas for the 
atom phase (cm > 0) one employs a — and 7^,(0,0) = 



For the classical propagator 



Vtj, = iq + 
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(C2) 



we recover the classical dispersion relation uj = q 2 /AM. 

In the broad resonance limit, the molecule propagator 
is, however, strongly dominated by fluctuation effects, 
generated by atom exchange processes. In the molecule 
phase (tM < 0) care is needed for the correct atom prop- 
agator Pp 1 17 . In a nonvanishing magnetic field the en- 
ergy of the atoms as compared to the zero level defined 
by the molecules is Ea = —&a + £73/2 + q 2 /2M. Here 
S accounts for the level splitting between the two atom 
levels (T3 is the Pauli matrix) and for a chemical poten- 
tial related to the difference in atom numbers in the two 
levels. For an equal mixture of atoms in the two levels 
one has 6 — 0. In the atom phase one has a a — 0. 

We can now write P<f,(iu,q) = -u+q 2 /AM+AP^+rh 2 ^. 
We concentrate first on the "molecule phase" em < 
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17 We neglect possible modifications of Pp due to a nontrivial wave 
function renormalization Zj, and gradient coefficient A^. 
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Here, the last piece subtracts the piece for q — 0, to — 
0, 5 = which is already included in the definition of 
rfv^. In particular, it contains the renormalization 
and therefore provides for an UV - regularization of the 
momentum integral. With uj — M(uj + 2a a), 8 = MS, 
qo = 2Mq' one finds the UV finite integral 

A ^ = -^J-0-J d q Q { [ {q > - q/2 f -Q + S + iq Q ] ^ 

x[(q' + q/2) 2 -w-S-iqof 1 
-[{q^-2Ma A f+ql}- 1 }. 

The qo integration yields 
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and we obtain for uj + 2a a < q 2 /(AM) 
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For uj — the result is nonanalytic in q 2 , AV$ ~ 
^fq 2 . For negative uj, however, the infrared divergence in 
dAP^/dq 2 is regulated. 

For uj — q 2 /(AM) both AV^ and the classical contri- 
bution to V$ vanish. This zero of P^ IjClll then defines 
the dispersion relation for the molecules 
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We recover the dispersion relation for fundamental point- 
like nonrelativistic molecules of mass 2M. In the broad 
resonance limit the composite molecules just behave as 
pointlike spin bosons with mass 2M, even if the mi- 
croscopic interaction between the atoms is completely 
pointlike and truly microscopic molecules play no role. 
(In this limit is approximated by AV^.) 
We can recover the relation 1C7|I by an expansion of 



AV<$> in q and uj where V<^ — 
A, - *■ 



-Z^uj + A^q with 
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The dispersion relation reads now uj 
coincides with eq. l|C7ll . 



(C9) 

(A^/Z^q 2 and 



The dispersion relation is independent of the binding 
energy cm = 2(7^. Hence we can use the dispersion rela- 
tion even in the vicinity of the resonance, cm — > 0, where 
both A^ and Z^ diverge. In our approximation the dis- 
persion relation is also independent of h$. We observe 
that in eq. (|C8(I we have defined Z^ by the term linear 
in a; in a Taylor expansion of V§ . The result agrees with 
eqs. (|3(JI7(jp> where is defined somewhat differently. 
On the other hand, for the atom scattering discussed in 
sect. [DJ eq. <|D1[1 we rather should employ a definition 



Z^Uj) = 



j=-cm ,q 2 =0 
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which differs by a factor two from eq. 1)70(1 . Of course, 
the dispersion relation (|C7|) is not affected by the precise 
definition of Z^. 

We next turn to the "atom phase" for cm > 0, where 
a a — and > 0. As long as uj < q 2 / (AM) we can em- 
ploy eq. (EH) with a A = 0. For uj = q 2 /(AM) the fluctu- 
ation correction A V^ v anishes. For uj > q 2 /(AM), how- 
ever, the integral (|C4|) becomes purely imaginary such 
that 
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Indeed, the correction AP§ accounts now for a nonzero 
decay width of the molecule - the resonant state can de- 
cay into two atoms if its energy uj exceeds twice the ki- 
netic energy of an unbound atom with momentum q/2, 
i.e. uj > 2 ■ (q/2) 2 /(2M) = q 2 /_(AM). At the threshold 
uj = q 2 /(AM) the decay width T vanishes. 

Obviously, it is not possible to find a value uj(q 2 ) where 
■p0 vanishes identically. For uj < q 2 /(AM) one has AV^ > 
and therefore P$ is strictly positive for > 0. For 
ui > q 2 /(AM) the nonzero T prevents a vanishing P^. 
The location of the resonance may be defined by the zero 
of the real part of P^. Similarly, we may define Z$ by 
the uj - dependence of the real part of P$, i.e. Z^(uj) — 
-Re(P (j) )/uj or Z^ = -dRe(P 4> )/duj(uj = 0). For the 
atom phase, our approximation yields Z^ — 1 since AV^ 
is purely imaginary. Associating the "binding energy" 
(-m > with the location of the resonance defined in this 
way yields 



(C12) 



We note that for a resonance the precise definition of the 
location of the resonance remains somewhat arbitrary. 
For example, our definition of cm does not coincide with 
the minimum of \P^\- Due to the uj - dependence of T 
this rather occurs for uj = max(o; m j n ,0) where ui m i n is 
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FIG. 13: Scattering with tree exchange of a molecule. 
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For q 2 = 0, this yields uj mm = e M - h^M 3 / (32ir 2 ). The 
maximum of the cross section is given by the maximum 
of T(s)/\P^,\ 2 (cf. eq. ljD14p ) which occurs for u> = s max 
given by 
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f 2 - e 2 M = 0. 
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APPENDIX D: 



ATOM SCATTERING IN 
VACUUM 



In this appendix we discuss the scattering of atoms in 
vacuum 18 and extract the scattering length. Once the 
effective action T is computed it contains directly the 
information on the one particle irreducible Green func- 
tions including corrections from fluctuations. Therefore 
the scattering amplitudes can be extracted without much 
further calculation. Contact to the observable parame- 
ters can directly be established in the T — 0, n — > limit 
of our approach. 

The amplitude for the elastic scattering of two atoms 
corresponds to the effective four-fermion amplitude. In 
presence of a field for the molecules we have two contri- 
butions: from the tree exchange of molecules and from 
the one-particle irreducible four-fermion vertex. The tree 
exchange contribution can be obtained formally by solv- 
ing the field equation for (f> for nonvanishing ipetp. This 



18 The "vacuum" means here the state with n — > 0, T — * and 
includes the presence of a homogeneous magnetic field B. 



solution becomes then a functional of ip. Reinserting this 
solution into T yields the "tree contributions" to the four- 
fermion amplitude. 

Near the resonance the dominant piece arises from the 
exchange of molecules, as depicted in fig. With q\ = 
-<f 2 , q\/2M = ql/2M = s/2, q = q x +q 2 = 0, u = uji+oj 2 
one obtains an effective interaction 



A e //(s) 



h 2 



+ K 



h 2 



Z$(s - iT(s) - e M ) 



+ A„ 



(Dl) 



This formula is valid both for positive and negative cm- 
The first piece involves the inverse Minkowski - propaga- 
tor P7 (iui, q) which is discussed in detail in app. El For 
positive cm the molecule state corresponds to a resonance 
with inverse propagator P^ = — Z ^uj + A^q 2 + rhi+ iF(s) 
and 0Ji t 2 = s/2, u> = s. In this case one has <|77[) 
= v = Z^cm and we use T(s) = T(s)/Z c j > . For em < 
the vacuum corresponds to a = <ja and the inverse 
molecule propagator becomes P$ = —Z^oj + A^q 2 -\-'\T{s). 
On the other hand for the atoms on shell (with a — a a) 
one has now wx,2 = (s — ejvf)/2, u> = s — em, resulting in 
cq. (IDTI) . 

For low momentum scattering one can neglect s. For 
the molecule phase with em < and for large we 
neglect the "direct part" from A^o and find 



M (O 2 



47T Z^CM 

Comparison with the definition of a in eq. (|16|) yields 
a v — em 



(D2) 



a Z^cm 



(D3) 



Here we have to evaluate Zq at u = —cm, 5 = 0. Using 
the result (fHTO)) in app. Oone finds 

Z*{u = -e M ) = 1 + ^ M 3/2 e M 1/2 . (D4) 

47T 

This can be compared to v — em extracted from eq. I|75|l 
_ _ hi 



v ~ e M _ "0,0 11^3/2-1/2 
- i7T M £ M • 



(D5) 



We conclude that for large /i 2 i)0 Af 3 / 2 e fl/ 1 ^ 2 the physical 
scattering length a coincides with a 



a = a(l + 47r^T 2 M- 3 / 2 e^ 2 ). 



(D6) 



In the broad resonance limit we therefore identify a with 
the physical scattering length. 

We find from eqs. {75J|, the familiar approx- 

imate form 



a(a A ) = a bg (l+ B _ B ) 



A 



B-BoJ' 
o^(M/3)-V>. 



(D7) 
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For 6 Li this yields A = 300G = 5.86eV 2 , with a bg = 
-1405a B = -0.38cV _1 . 

We next turn to the atom phase where cm > 0. For 
a first discussion we also omit the decay width T which 
vanishes oc y/s for s — > 0. The effective interaction be- 
comes pointlike 



A 



eff 



Z^tm 



+ Aw, o 



(D8) 



and we see that Z^cm replaces (D\ — 2fj) in the "mi- 
croscopic vertex" Q. In the atom phase the physical 
scattering length reads 



a(B) 



\ eff M 



An 



' AnZ^ M (B) T An 



(D9) 



The first part corresponds to a res (recall Z^em — v — 
m^) and it is suggestive that the second may now corre- 
spond to a^g. 

The situation is less satisfactory as for the molecule 
phase. In fact, a = — h 2 ^ M/(AnP) accounts now only for 
the resonant part of the scattering amplitude, resulting 
in the relation 



a=(Mpy 1 / 2 (B Q -BY 



(D10) 



In the broad resonance limit the interaction is pointlike 
and we may absorb the background scattering A^.a into 
a more complicated B - dependence of h^. This would 
result in the relation l|D7|) also for a in the atom phase. 
The B - dependence of a or c therefore depends on the 
splitting of the microscopic interaction into the two pieces 
A^a and —hi/pA. Since all our results depend only on 
c the B - dependence of the observables is now affected 
by an unphysical "Fierz ambiguity" 60] . This is clearly 
a shortcoming of the approximation, in particular of our 
neglection of the molecule fluctuations for the computa- 
tion of Z$ and the renormalization effects for h$. The 
uncertainty arising from the Fierz ambiguity is reflected 
by the two curves in fig. |3 where the solid line reflects 
the choice l|D7|) whereas the long-dashed line corresponds 
to eq. (|D10|I . Presumably the solid line reflects better 
the true answer since at the end all macroscopic results 
can only depend on the combination a res + a^g. Close 
to the resonance the Fierz ambiguity becomes irrelevant 
since \a,b g \ can be neglected as compared to |a res |. 

Close to the resonance the scattering at zero momen- 
tum provides the same information on the relation be- 
tween c and B as the binding energy 6m(B) in the vicinity 
of Bq. Indeed, the dependence of the resonant scatter- 
ing length a res or the cross section on the magnetic field 
involves the same combination of parameters. The res- 
onant scattering length a res = (M/3) _1 / 2 (i?o — B)~ x in- 
volves just as the binding energy in the "molecule 
phase" for cm < 0. From the experimental parameteri- 
zation of the scattering length 



a(B) 



a bg 



1 



A 



B-Bc 



(Dll) 



we can extract the parameter tb as 

tb = -(abgAkFy 1 - 



(D12) 



This does not involve the Yukawa coupling. Experimen- 
tal comparison of cm {B) and a(B) on the respective sides 
of the Feshbach resonance therefore amounts to a test of 
the computation of the vacuum properties. 

For a realistic scattering of atoms we have to take into 
account the momentum dependence of the atoms and the 
nonvanishing decay width f . From app. |Owe infer (for 



and to 



e M 9(-eM)) 



f(s) 



^2 M 3/2 s l/2 
47T 



(D13) 



On the atom side we find a resonant scattering for s = em 
with a Breit-Wigner form. 
With Za, = 1 and 



Im(A e// ) 



(s-e M y + f 2 (s) 



(D14) 



we may compute the total cross section <r tot as a function 
of s. Using the optical theorem otot = M -1 I m (A e //) one 
obtains (|u| = \q\/M = y/s/M) 



O 'tot 



-Im(A e// (s)) 



(D15) 



4tt(Ms)- 



1 



167r 2 (s-e M ) s 



For s — > 0, €m 7^ this yields the expected result 
M 2 ht 



a tot 



4ne 2 M 



Ana 



(D16) 



Independent experimental information on Jim can be 
gained from scattering in a momentum range where the 
pointlike limit does not apply anymore. In particular, 
the peak height of the resonance at s = cm depends 19 
on JiA, 



An 
Me 



M 



An 
AID 



167t 2 or 
M 2 h 2 



16TT 2 h~ 2 M- 5/2 p- 1/2 (B 



B ) 



(D17) 



The value extracted from resonant scattering i|D17|) 
will again be a renormalized coupling in vacuum, but 
now evaluated for the particular resonance frequency. 
It is therefore conceivable that it differs from /i^o, eq. 
We note that the cross section remains finite for 



19 More accurately, the contribution from the background scatter- 
ing a^g should be added. 
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all | cm | 7^ (B 7^ Bo). On the other hand, one finds for 
cm = and s^Oa divergence a to t • 

It is interesting to look at the behavior exactly on the 
Feshbach resonance (cm — 0, B — Bo) both from the 
picture in the atom and the molecule phase. With 

h 2 

Kf f (s) = -pj^> = - z <p( s ) s + if (° 18 ) 

one has 

|7^)| 2 = Z 2 ( S ) S 2 +f 2 . (D19) 
In the molecule phase T(s) vanishes and 

Z ( S )s = s + ^M 3 /V/ 2 . (D20) 



On the other hand, the atom phase is characterized by 
Z tj> = l and f (s) given by eq. 0I3)l . While |^(s)| is the 
same, the phase of A e //(s) jumps between the molecule 
and the atom phase. This is a typical shortcoming of 
our approximation where the molecule fluctuations are 
neglected in the computation of Z$ and T. In particular, 
the molecule fluctuations will provide for an imaginary 
part of X e ff(s) also for the molecule phase. 

Our functional integral approach reproduces several 
results obtained previously by quantum mechanical two 
channel computations in the vicinity of the Feshbach res- 
onance. The fluctuation correction to the mean field 
binding energy coincides exactly with the result from a 
quantum mechanical coupled two-channel approach |6lj| . 
The estimate for the Yukawa coupling given in |35| coin- 
cides with our "vacuum coupling" h^o (|96|) . 
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